{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Contents" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook covers the basics of creating TransferFunction object, obtaining time and energy resolved responses, plotting them and using IO methods available. Finally, artificial responses are introduced which provide a way for quick testing." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Setup" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set up some useful libraries." ] }, { "cell_type": "code", "execution_count": 39, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "from matplotlib import pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Import relevant stingray libraries." ] }, { "cell_type": "code", "execution_count": 40, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from stingray.simulator.transfer import TransferFunction\n", "from stingray.simulator.transfer import simple_ir, relativistic_ir" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Creating TransferFunction" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A transfer function can be initialized by passing a 2-d array containing time across the first dimension and energy across the second. For example, if the 2-d array is defined by `arr`, then `arr[1][5]` defines a time of 5 units and energy of 1 unit.\n", "\n", "For the purpose of this tutorial, we have stored a 2-d array in a text file named `intensity.txt`. The script to generate this file is explained in `Data Preparation` notebook." ] }, { "cell_type": "code", "execution_count": 41, "metadata": { "collapsed": false }, "outputs": [], "source": [ "response = np.loadtxt('intensity.txt')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Initialize transfer function by passing the array defined above." ] }, { "cell_type": "code", "execution_count": 42, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(524, 744)" ] }, "execution_count": 42, "metadata": {}, "output_type": "execute_result" } ], "source": [ "transfer = TransferFunction(response)\n", "transfer.data.shape" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "By default, time and energy spacing across both axes are set to 1. However, they can be changed by supplying additional parameters `dt` and `de`." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Obtaining Time-Resolved Response" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The 2-d transfer function can be converted into a time-resolved/energy-averaged response." ] }, { "cell_type": "code", "execution_count": 43, "metadata": { "collapsed": false }, "outputs": [], "source": [ "transfer.time_response()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This sets `time` parameter which can be accessed by `transfer.time`" ] }, { "cell_type": "code", "execution_count": 44, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([ 0., 0., 0., 0., 0., 0., 0., 0., 0.])" ] }, "execution_count": 44, "metadata": {}, "output_type": "execute_result" } ], "source": [ "transfer.time[1:10]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Additionally, energy interval over which to average, can be specified by specifying `e0` and `e1` parameters." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Obtaining Energy-Resolved Response" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Energy-resolved/time-averaged response can be also be formed from 2-d transfer function." ] }, { "cell_type": "code", "execution_count": 45, "metadata": { "collapsed": true }, "outputs": [], "source": [ "transfer.energy_response()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This sets `energy` parameter which can be accessed by `transfer.energy`" ] }, { "cell_type": "code", "execution_count": 46, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([ 0., 0., 0., 0., 0., 0., 0., 0., 0.])" ] }, "execution_count": 46, "metadata": {}, "output_type": "execute_result" } ], "source": [ "transfer.energy[1:10]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Plotting Responses" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "TransferFunction() creates plots of `time-resolved`, `energy-resolved` and `2-d responses`. These plots can be saved by setting `save` parameter. " ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD5CAYAAAA+0W6bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvXeYJVd55/85VTd3DjPTkzSjNJJGSEISiCAUjAUWwYRd\nk9bYYLC9XgdwxGDv2rLXYJDX4Ih/uyaYh7iYZFhEELYGJEAiCUkohxlpRqOZ6Znp6Xxjnd8f7zlV\n1a0Ot3u6+94evZ/nqefeW7fCuaHe89b3fc97QFEURVEURVEURVEURVEURVEURVEURVEURVEURVEU\nRVEURZlNIcBC08vxtWybWcuTKYqirEPsdedkmtrwuvvrsIZ2NVirEynrhn8B/merG7EA/w04DIwB\nfS1uy0pwGjCOOlNtjTHNLWuNGvBTixzwQWAfYuBuB65d4jH8reBs/ggxNOPANFBPvb5rec1dMlng\nr4GfBrqBkRU45h7k84ynlmetwHHnYx/w/NTrx4Au5v7OlTYhaHJpRbuUU4cMYhCuRAzcfwc+DexY\n4nHm8iXehRiaLuDXgO+kXl8wa9/V8kWGgAJw7zL3n+v/boHfIPksXcBtyzx+M1jU2153qAeurAVT\nwJ8hRhzgy8Be4JIF9rkY+BHisX8KMZCLMdtI7wH+Avg2MAmcAfwScI877sPAr6a2vxo4APwuIocc\nBN6Yev/FwN1u3wPA7wFnA/e5908A33DPzwVuBI6591+VOs6/AP8E3ABMuPM2yx7gzanXbwRuTr2O\ngP8KPIDcCfzDrP1/heTz3418zx9FJJMvIZ7+7wM73bH8tbgF+KL7PA8Cv5w65nVIh/wRd9yfAJcu\n4TMpy8Q0uaw1asBPbTYBuxADMhc54AuIQegD/hX4zyzvdv71iLHpBB5FDPNLkDuBXwLehxixdNu6\nEYP1ZuAfgR733gcRg98NnA/8B2LMdrv3e4BrgA7EeH8M2AC8Fng/cF7qPK9DNP1OpIOZi7muvfmk\npDQvAZ4BXAi8GvgZt/5VwJ8Cv+A+w8sQg/wLSOf6UsTT/19zHPNTbpvNwM8hdz4/lXr/Z4FPIt/B\nF3lyx6GsAqFpbllr1ICfumSBjyNe6APzbPNsRHb5W6ABfBb4/jLOZd157kW8yTri9e51738L+Dpw\nRWqfGvDn7rxfQTzkc9x7VcRwdwOjiJYPTza0L3Xn+Ig774+BzzHTC/8C8F33vDJH2w3wd4gXPQL8\nYNFPm/BuxBPeD9wEXOTW/zLwHuCH7vXDJHdFC7EdeC7wh8h3cAfwAeAXU9vcDHwV+c4/ljqnsoqo\nhKKsJQFyu14GfjO1/iskgbr/gni/j8/a91GWdze4f9brFwG3Ip7nCCKLDKTeP4YYXc8U4iWD3AW8\nGAn47UE6mrnYgQQcR1LLf0G8exAjN7tds7HAbyF3IH2IR90sh+Zp/zbEaC+VLUge8WRq3WPA1tTr\nw7POWUCv41XnJCSU7Ujnfjcieb1l1vu/h1wH/al170DuOO8DXrhQu5pLblTWEwaRIDYgRrCReu9F\ns7a9ipnGAcQoPrSM86blhjzizb8e+DfXhs/TfMfwA+AVQIgY108j2vFsHgO+ySJ/8mUyiUg0nqEl\n7LsfOGue9xaSZQ4iF3InckcC8rkPLOHcyipwEt51Dfgd5O6wE7kruxG5W90OvABxmjy7gde4x61I\nrGcXM52dGO25Tz3+CQnsvYy5JYM030Hkjrcgkst/Ap65zPOm/+I5txxF/ngvonkjmwV+HtF4G8jd\nQmOebf8f8ud+vdsvi7T/3Dna1GzbPT9Gvo8iYozfPMc2s4/hj/MBJEB5iVt3FkkHdBg4c55j7Ed+\nk79EOsELgTchUonSQk7CAz+E/JdAOuV7kTstgPcCb5u1/cuRGEcNuQN9CLhsvnapAT+12IEE/y5C\n/jheLnndPNvXECP1RkTSeDXiOS/GXAG+9OtxpFP4NCIJvA7xxOfbfjavR7TtUeTz/Pw8+00gHcNr\nESnoCcT45RZo51zMtc37EB36MPBhxIjaBfZJn+szwDuBTyAa+edIBh39JZLeOYJk4cw+1uuQzJSD\nbr8/QYK4830ezR9fAwLT3LIIO5FA/m2IoT4A3Dlrmy3MvOM6wJPvkmM0H1VRFGVh7PsumFttfmgi\n4qHJpA/92hELc9vVTpJ0268juvgLkA5+LxJ7OQb8PRI7+rjb7wNIQsDn5jq/auCKoiiLMJ93vasr\nYFdX8vprR+pzbZZF7mw/hmRFXYB443e497ch2vizkDvJ7al9t/HkRIOkXc00XlEU5alMYGxTyxz4\npIJ7gL9x6+5CMqVOd8sBJF5yGMntfy0iA56ODGD73nztUg9cURRlEU5Ca74ciencSTKe4Y+QlF5P\n2vLfg8SO7kESDH6dBeIcqoEriqIsjP3Hi8KmNvyNOxqg5WTn56ILzl9KcXVddNHlKbxc9bznWlaA\ndh2Jue4klDvuuhs7OdzqZsRc987rue6PZ6dytg5tz8JoexbmVGuP6diwIu1oV6li3RlwRVGUtaaJ\nHO+WoAZcURRlEdrVgK+FBr6PJALr02H6kXoADyBJ7b2p7Zsu5NIOXH3F5a1uwgy0PQuj7VkYbc/c\ntGs98LU4516k6Hx6tubrkToZ1yOlM/uAtyMFXD6B1LOYr5CLbScNXFGU9sVp4Cdr5+yHLm0uC+VN\nPzw1s1Bmf6CXITWccY+vcM+XVMhFURRlLWhXD3wtDLhFPOkfINNMgYxC8nWND5PUb15SIRdFUZS1\n4KmcRng5UiVuA6J73zfrfZ+zOR8rksepKIqyXNo0hrkmBvwJ9ziMFPW/DPG6h5CSp5uBI26bpgq5\nXPfO6+PnV19xOVdf2R6BDkVRWsueb32bPTfPN/Xp8mmFd90Mq92sEjKryjgyu8nXkVnTr0FKJ74H\nCV72MjOIeRlJEPMsZnrhGsRUFKUpViqI+fHLmgti/vz31jaIudoe+CbE6/bn+jhixH+AFGx5MxKs\nfLXbZkmFXBRFUdaCNnXA27ZdC6EeuKIoTbFSHvinntWcB/7a204tD1xRFGXd066erhpwRVGURWjX\nofRqwBVFURahXetuqwFXFEVZhHZNI2zXjkVRFKVtCE1zyxxsR2agvxv4CfAWt/5Vbl0DmQ8zTdMF\n/dQDVxRFWYR5JixuhhrwO8CPgU5k9vkbkYmNXwn871nb7wZe4x7nK+iXtGu5rVIURXmqEDS5zMEh\nxHgDTAD3IjWf7kPKac9mSQX91IAriqIswgoVs9oJXAzctsA2SyropxKKoijKIqyAp9sJfAZ4K+KJ\nL4V59Rs14IqiKIswn3d9xwnLnaOL6uNZ4LPAx4AvLLJtUwX9PGrAFUVRFiEzjwG/tM9waV/y5sce\ne1Ks0QAfROo7/c08h08f/YtIQb/3ItLJ2SRTUT65XQu2WlEURTmZPPDLgdeTzAsM8EdAHvh7YBD4\nsnvvRSyxoF+bpqcviBazUhSlKVaqmNWeq5pTwa/+ZrQS52sa9cCVVcFOj0BtCjs5DPUyNCowfhh6\ntoKNIMhAkIWoLs/z3RDmMJkC5DogU8CE2VZ/DEUB2nckphrwUxxbm8Z+/FexkycAqI2e4GM3PMih\n6WSbjizkA6hFEAYSce/KQmShmLEM5KAaQW/O0rAw1GGo1S3bdm3iyN5hejd2Mj4yzYYztxBks4SF\nIo1KmVzfICZXoDE5Rti7ETtxHNPZL1dDLg/VijSgUIJMDoIQawKIGhA1sGEGsgXId0AYQiYPo4eh\n0OWOUYTyhDxm8nKMqREo9srr2jSEWSiPQf9pMD0m56tNgwnkvUKnnK88Ic+DUDqUqA6FbtmuMgGV\nSQgz7vWk7Bs15HVUl8dsHmoV6aBqZch3Qr0ibQ1z8hnqNdk2k5XHIJTjhFn5rNYmnVqYlf2D0LU3\nl6wHeS9blPdA3vfrM3kI89CoQlQDEzor5NvijmEjd76s7IORztbaxGqZQNZj5THIYHybrJVjhNnk\nuX/PGIgi2S/MJu00IdiGHMtGblvv4fpzGjnObNLbBSHGrE0mdLvmW6sBPwWwUQOmj2OP3s/fvPrn\nqDSgEhmMk87yIeQC6MhAZxaKGdjVY+nKQF/O0p+HRyYMl5zZQSYbkslliBoRWEv3eRfSmBynMT1F\nbtuZADSOHyI840KoTHH6Mwz0bKJv9DBs2AljwxBmyGw8SwzW8F5CZySNjWCDHINMHqqTMHFMLvz+\n7TMNWq4kj4H7i4Y5qE7A5t3QtRnGD4qBKvQ4Q1aC6eNQ6BWjYBtivMK8M3whsXEK8+79mhw7yEh7\nTOCMUR7KI3LMICv7mBCTLWEbVWlHmJX93N2EyXZga1OJ8Q1C+XzZkuyfKWCCLLY2CYDJdSYGUH7F\nuCOwtSkxckFGDJQ3zNjE4NlIPk/agAXuedRwK4w7pnGdUsMZVfdonIE0KUMahEBn8t3ERtedM2ok\nv5ExrlNzx/SGGus6ooY8+u++XnHrIteJZPwfWNrpjx9mxPDPSK42ybbx5w3dP9zI72ntjPdWUslo\n12qEbdqsBVEN3BHtv5XPvOHl3HvCYN1PGRpLdxYK7trszEJHRox1zcKWoqWjELDlgjMZP/AEXVuH\nCAtFrLVkN2yFfFG8xHwx8XB7NopB9d5urSwGtzoJHQNy0U2fkMeBHTMbOT0GHX3iKQJkConXOnVc\n9s93J4baG5/yCSgNwPSI7JspJB5lmBOjmOuE8qhILt4YRTVMoQ9bHpH3a1NgAjGw9emU9+cI85gg\nxEZ16WRcO22jir88TJBBjFMOrMXWpmS7tPG1VjzoWllWBWEiFXkPu16R12nP1nus3rAFoRgvGyWe\nuzFiUL2RzRTEq27UpI2ZvGtGIzG41s0VbiPphOL3UoYb3LHdnVCQS4xxVEt+EyyYjKzzx4sNu9ve\n1p2Xn+oU3O8h57eugzGJoU4fY7ZH779T30n6uwf/HcSf1//u9Zmfyf02Qd9O4h9y+djv/FRzPvhz\nb1INXFmE6Ja/47Z3/wU3HgyI3M1daCwdGTHc3TnxrMsN6MnCUMGycUsn0+MV+rf0cezxETKd3fSe\n3UGQyRL0D4lxLnWJDBFWEzmha1CMT60MG06HsSMwuFOM0eDpsj7MQs+QGJJ6RSSKriGYPArdm8QA\nN6rJhVfoEV3ce3xBRvYDyHfJxV3oFe+1NpV41QDFfjHu1opBz3U4D9V5YWEOWx2HbEdiJIIMNqpj\ncl3YRgVjctj6NCZTdAa2Jp5u1h0nqmFMSkrxxth1HsZLEN7oGv9+XYx4w+3jJYdMRoytN4jWyneV\nNt6QyBBhNjFCtgGEMw1/VBeDmnWeua07gzVLkvDSCYhxjretJ222uN8lcOuzKWPqjKLJJJ2Jl1bi\nuwDkP2Od5x2kPHDsTE84qss677lHkRh1L5vZWtKB+c+f7rjSQrQ/d+Q62iB9N0Py260Q8xSqajlq\nwNcR0dH7+duXXsVo1QAhobH0ZS0FJ292ZOCMzohyw9CRsYQGBnuybH3WxdTHx+isVshv3EzHWTnR\na7sHqD9yF8HGHXDkUSh2y4HOejYcexTGjsLQ2VAeh+6NzhsfkousbxtMDEPXpuTiqkyI19yoycWT\n7xRDmik4XTWS1wD5HrnIigNiaDBi6Cvj4pHnOp0XviGRLkyAMQab65R9ohom34OdPALZIiZbig2W\nbVQxYS4x3JVRiEIx2o0qJtvhDKqV9tUriVHId83sUHzn442RN9yZQur234ohazgjFISJnu31dmmY\nHCPtpPljZwqJ7GAj8AWUvIRio5RH69pq5L8Qyycg5/QSjyfuKMLkeL7d4OQMk3i1NhJjHnv/NnU+\nZp7fd0BB6Lz0lCGNDX5K247Pmf5O3PdgQpHejIF6NZF/0sdMyz/ehNlGopl7Tz5YOfPWrhKKGvB1\ngJ0+zo2vOI/vDMstqEGMc28e+nJQj6A7axkqWmoRXLK7l4mRafpP30r52FGqR49gMhlyAxswHb3Y\nyRMShIoaZC64Uk6ycQeMH5OL5vh+2HiWGOAwJ7JIJg+dGyQgWOgWr7g0AOOHnIddg+4tMDksF2a+\nB6rjzijVE23Za9PZkiyTR8Srzpbk4s91iLEPXNDSNpxHXcLkOqA2jcl3YaePgwmxUQ3TsVEkkJQm\nbHIdsSRiyyOYjo1yZxB7gM7TrFeBKLkTMIHcVWRy8l5tOjHMOK86CBNZAlzHVEqMiDfqGGZ48UEO\nbJjSrSMgcN5jJOfzRsg23Oe3yTlsAzJF9xmyiYSS6xA5y3cwmULq+HaW8Q1S+rbzXKOUtm2t6yQC\nOV+UClJ63TvW0/3ncN+Dl4hA/jdRLfUvNkknEzgZyzaSTCR/1+A7EK+FBxnXgZrUNu7z+P2CEKy/\n+/Aef5SKBZw8J1GNcFVRA97mRLf8He9667uo24BcIN52Tw6ygXgFOzsiqpHh7KE8R4+X2TFUpFap\n072xBxNm6Nh2GmFHFyZfEq87V8D0bpTn9VoilUwcgw07xBj3DMHUqDxCEuTL5OX96ROyz/QIFPsA\nC41ALsYgIwalXoae02DqmBzDpQnGt/pe9igNYDJFCe7lOlPebka869qk6NQmcMbWOOMdiJRhJFBn\nch3OU2xgGzVMvjM2hpILPEua8IYgk0tJIe5OIpMX45DJidGtTcqdQ3065Yl7/GdxkkHUSDqtwHmY\nJnDvAzjZIszK9x96zzgdwHOGDee5ei0dHzCMUp2FTe4WfMaH9+7DbPJ9xlJPPZEb4uwVpzN7+aZe\nTTznuP2prJB0IDF+7jqOdKcWyx42OY7viHznhp25ve8c/P9uRrB2VszBy0q+3V5y8t/XCmaotKkD\n3rbZMQpgJ47w7t9+F3UrXnc2gFIGGlZS/kBkk60lizGwcUORnq0bKfb1kO3tJ9vTS6Z/EyabTy7Y\nrgHRaHMlSWUDMdQ7LhHjXOx2J49gehw6BsUTDzISdGzUoHOjaNNY8Z5zndC5SV53bHQGvyAG3l+U\n3hOtTWI6N824yKyNMLkusFYkDhOIHBJkMPluaFSwcbpaDlMacMG9AFudSCQFF8w06Ys3Dox5A+5v\nzyvOwGRSGSpVMR4+UJoO/nmDGR/ba9SpoFps5LzMUE9u9+M0PhBjXU+kjDjNMOd09GKiDXsD5Y1a\n2ntu1JL1ftvYqEVP1t29tJJO1fMBxch9Dr+NlzPi7JN6qiOYtS5tYH3sAJt8b+nAqv89/B1K/LsE\nyT7+db0i2/iOF2be2XjSWTDpTmEFCUxzy1qjBrxNsWMH+evnX0g1MhRDMd49OcnF3lqyXL4h4hXP\n28zOM3qJLGTzGXKFDOXjI+QGN5E743xM7ybYfCa1kWHYuksOPD0OZ14mMsHg6ZK3XJ6A6pRIJaU+\nuThr02LYS4MwdkiyPUr9Ip1UxuUCy5bE2FUnZzY+1ynasd8mzMmtv0uJs5XxmR4XYKNafKtssh3Y\n+JjGpeBlsLaBrZedNFKSd4v9iRfoDEScQZLWlb2x89p1Jk+SgWEk+Bhm3XsFMUT1cpKHjUk8ap8l\nYcK4I5mRtgjOcGeS5z5gmc7iCEKXqhgmhi/MuM5R7ixo+EBcyouFpFPMermEVMZGIHcP3iB6+QES\nOcuYJNgahEmHMcPXTBnXMJto4dYmOekzPlvaa3ffmdeqfVDTG98wmwQZTTizE40DoOkURZP8BnGA\nmJnetu/oguyKyickZ190WWtUQmlDom+8i0+/82+ZqAf05Ww8qGYgb7lkS4ZHj9UIQ8OxR48wPl5j\n5/lbsPU6ma5u8lt2Ek2MEh15jGBwK5w4THbnbjnw4DZ5rFWg0wUbe7fIusqkaNv1ilxcWy+U9ROH\nxagXehIvKVuUCzjfLcbb52I3quKB24as91qsD26FeZe2lxFpxMkfSV6wk1UyObkYjMFWJ11qn5FR\nmiDtDnOSURIPcsmKxBGEmEJPEqAzjSRAV5tO5B0/aCU9iMSnGDZq0uF4A+KNuU/hazg5Iqq5QJvb\nP9fhDFxdtquMOwPtUvDqLgc9bWhtI5ENGtXEYHpN2pDSp1PShc/vTudrx0FJJ9Nk8k6rdgYv8kbR\nHSMwiRQRH9N7y5F0QD6oGncEYeKl28jJNul9nWTidX8skPqOAZ8tlBhZZ3jDHDQayXeSzmOPBwa5\nu4Ugk9wVRfWkE/ASzQoP8NEgptIUtjbNX7zt74gI6M6K8c6563h7hyUIA56+ux+CENuoky9WiGo1\ncn0DZLq6sdPjBMUOKHbJSMcN28VgRw2X7mUlW6RekcySYm9yQZbH3YhA5y1PHIGNu2HscTFgxX65\ncEobZNBMdUK87ckj0LNDBr9UJ1JerPOYyieg2CeGu17G1muYYj+2Mi7ZIMZga9MuH7uGaQTY6WOY\n0gZMpoCtTsjAF6c323oFE2YxtpEEzoIs2LJo1kRggyRfPK1TB5nUiMVyEnz0erM33lEt6WACFyhM\nB9Pi/GqT6N4+AGgyyXcTSw0pKaM+LefwHmu9Kp1MOuCbfu49dLzxSud6O3yAMnRBP+8R2yg5xuxM\nGP85alOuY6vIObPFxEufocO7/dM6dZgjHrzkZQxIjKmPe/i7Fm/o/SAtL/XEv0822W52Dns8iCjl\noftH76V7rT/dGa0A7ZpGqBJKm/G5a3cSuUyTDnfNZ4ylK2vZeVoXX3+4xme/d4K77h4m09lNR38X\npW07qRw9TJArUJ8Yh4IYRUpdcGgv7LhYDrRlt6wvdicXc+eQG3rtZIUwK4HJbEkurNq0M+ipW+KJ\nQ1AZg45NYgTzPTA1nEgq1g2IqZcT3TpTAAwmW3KyRzJIw9Yrsn2uM+Xt+RzxEFPoTd3GV102yhRk\nCljbgGwBO3WMeECI11i9ph3VE5mnOpkYlUxe2ghODvHDz8vu7sEFZutuaDmBHNN7ymEmkRNiT9Zr\n0c5DbVSSc3ivNfAZJM4QxWlxqSCdMdIZ+DuFGQbJJJJLHIz0unc10bnTBtQfo15J9k2PckwPkIq9\n4VpinMEZfptqd0oTB2akCPrgY1qzhyQQ6j+fCZJRsD4YCzPvDNLGOx3UbFSJg6/xCNNUfvoKFjBp\nVw1cPfA2IvryH3P3iKEnZ8kHIpmc223JZw1hRnKgX3ZhB52bN9GYnqKwaQuN6Ukqw09QGNpKbeQo\n2V2XQEevDGm3Fs68BA78BHo3izEu9YoXlytB3044eIdkpPiLGyAzDljJ8Y7qSU0RLFTGMD3bseUT\ncrE7D1q8+LwUyQRsZQxTGhQ9OtcpHnamgK2OizGPg43uQo4y2PqYGOd6GVPokwu27kf5OY3TG6ds\nCWwUZ5+YTH6W1+pusW2UjFQ0zsutThIPXsm4EaL16UTC8ZIJyHa5DuL8Zt+R+EEw2ZJ0Zpmi+/5S\n2zScbOTbHdWdRFOX7602lWjxxkA9nb9tkrYFoWQE+hon8QAf12F5Kcan+tWrrq5KkEgQcQDTY2el\nE2blriFbSqXnZVIeswu6eonNG1W/byyZRMn3BklnGq9zEou1SSzA/2b+bsgyMxhsAjlu3FmlMmbi\nlMKU8Q4yxOmPK0SbOuBr4oGHSK3bL7nX/ciszA8AXwd6U9u+A3gQmfDzhWvQtrbBWsvn/voDBAaK\nIfTlLb056OrK0tFTYMOOQaYnqgy84pfIP++V5PoHKR8+SFStUthyGpnObkwYJkO5swWRTLacLxfq\n4YdkOe1SuQXv3gQnHoPBM6FvhxiUriFZ8l2JVpvJi0eeKcjrjg1SqwNg6qgEE61NZAenU5t8t0gf\nac/SGBlU4/KKbVVmlhIZxWWe1MvO6KXS6tJeVdRIOpo4t9x5mJVxMQC+JgkmGdnph4JPHU887SAL\nUTXlqUby2mv93qCnJRATEnvi4IxqUTqZ2JimvMTYI20kAbZMQXLk4wE61nUCTo6oTycDi2rTiXGM\ndWGfZuiDhu4cDefRzk4T9AFPH0BMp97F79kk1dMbQ58v7z+3zwf337eXdvx34UdOxsbUJt/l7ACp\nl3bSwVnvUccjKNPG2adumqRTSEskqYD4nJkqJ4kxpqllrVkLA/5WpDi5vz98O2LAdwH/7l4D7AZe\n4x6vBd6/Ru1rC+yX/pD7Rg1DRTHgWQODeUvf5l6GnvNcTEYKTE1+41Mc+cQ/su87d2LrNcJiCdO3\nierxYTJDO6kfeAgmT8DW82H380XX3rIbLnwRnH+tXGBDT5NBOIXuxOOqTsigHK8LB6ELSLqfbXJY\nLpLyKER1TL4H07tTBtDUy+JVe2PqtdIwL4NwMJhMXgy2kWCkLZ/AdA5BPAQ+I+cOsmJofWCrVia+\nGL1R83VFvPce+OBqQfYv9BLXEKmMJ0FEb6iDUAy1r6cSF3NK5TmHeXdsZ6AarvKfbSTbe625NikG\nvl4mLt6UHmTiL+xcKTGA/rz+nJlioiuns12Ml4QaybF8xou1yfuQ8vKdxBFLPT4tMCXBpNMC486x\nPrP9fl067923Md6/kXw+fz4/+CaWS7yRn5V26QtQpTNjvMFPB09j6SSlpfs7gLSk44OZ6d9xhfA3\nSYsta81qn3Ib8C/AO4HfBX4W8a6vAg4DQ8Ae4FzE+46A97h9vwpcB9w665inZDGr9zxTBs0UQthc\nspzdZdnQXyDMBHRvHqT49OdJtT8gDvx0DFC76WNk+zcm3uqZl8CRfTA5CptOl9vp7kHYdlFSRjVb\ngNFD0L/DZUuMyUXnC051bHSGuiaGZnpELqKCBDxNaQN2ajjxcFwQ0QQZMdjeG3cZJ3EwzN86+yyD\neOCIN9x1rG3Eoyhjo5WuOxJXt3Pn9VqyDxg2KomxShd0Sq+Lsz+8Tu2MTVw9L0gMuHE6eMNVPPRp\nfd44pL1un2nj2+G3t/VEi06POgTiwUTpsrHpFDjfHq8JR/WUHOK+03jou0kCgiZIRpF6jXu2R+6P\nFWd0ZBL5aHawNq4JYxNLFd9pmJnP0wFKT9xJzsoJb1Rnesr+vxJv6/0+fx6bfFZI9PpYpkn91zAE\n3VvczieFffAl4eJbAWd/ubES52ua1dbA3wf8AdCdWrcJMd64x03u+RZmGusDyJxwpzzRx3+Z6YZh\nqChZJ88cCslkAzbsPouwWCLo3UT1vh8y+rV/o2Ogi2xPP/XJCYpnnEf2witkNOX4MfG8K1Ow6Qwx\nHMVu6OjyVsaoAAAgAElEQVSHww/C8MOSOljqg+lRl0ZYFUPde5o8es9p+rhcEIVeMcZ+Xb4nMcCx\n9+c9Ml8DJBlKbfyQaQIwocseSckSQVYq++W7oFGXAT2hBO1svYzJdWDLoxL0hNjA2coEpmPQedfO\n8Lr88Ph1QKoDKGCnjkp6odeDvRH12m6QAeMNp9dxM+KpB9nEq4szU9x3ZdPGLEgZEJN0OqELlnrD\nSsrgQapTscml74OK6ZGjvqOZK7AZnz+le/vOID4vKZ28nDq2Ec24UYMgSgUto+S8aWloRgfTSLb3\nHUgsm/g2puSPuCP3302qI7KR/G8z+cRux9+TdTZ51nfm/2Oxcbfy+5nmDG6znIR3/SHgJcAR4AK3\n7jLgH5DJjuvItGnfd++9A3gTUp/gLYjMPC+rKVG8FGn07czfI/n7n/lY6L1TgmhkHxO330I+kKyT\nS/ojNu3axqaLdhOVpznw3R9R2XcvuW1n0rlpgLDUSX1ijOOPHebxm26ifvd3KH/rc1TuvEWCk6c9\nXYx4sRuOPCLedueAaN49W2HyOGw4V6SJqC71S8YeT4Zdg3jk+W65eCYOiSdZ6BPdNgixlbG4Tkk6\nc8DWpmJv1ZhABt0A2Ajri1kZI8+dxmp8YMxGcSlXWS85zKbQk3xZrjaJCULRsr2X6dMIsyXA6fF+\nWL8JIKphir1J8SlvWGpTyfNYEvAyQoakRKrz2P0Qc0j0/vhzBXIH42UUE7o64PnEWPqJEvxdhU8H\nTA+Q8XcBvmRsHLjFGdM6cWaNlyV8emM8oMnlZ/taIvGdi79jIdXRBUnRKK9/+w4Aks+brgEee8Cu\nzd5ox4HFKDlPHHD1f/h6oml7o9uoJtvPlsfi/W3yO8XaeD2pVpgOlPo4if8/rwBBEDS1zMGHEUk4\nzfXA/wAuBv7EvYZlyMir6YE/F3gZ8GKggHjhHyWRTg4BmxEjD/A4sD21/za37klc987r4+dXX3E5\nV195+Qo3fQ350nW89z9G2VJytbo7c4w+foTaIwfZ9gu/xmnnHpLb4COPEuRyNKanqIxNEIYB1loy\n3X1knv7TcqzpMbh3j2SV2Ah6h6T8a+eAXBD+lnrsgOix2RJgZCCOz4/2+d5RA1Psw44/IUHMMCsy\nha8THeZF+/UFpjo3YX0ecaYogUqXemgbVSkq1ZgADKbQg61OYVLFp8Bixw5g+s+WWifTxzDd25Lv\nyY+qjPO68+Idp2+ffd3peIBQQzxknyNdnUhkAmMTY+SzVxpebw/jNEVZX4k7AvFmvYZeSxl460ZQ\nBonRCfOyLpYIDPhh9PEIwpRHbwxxnZOIJMjpUx9hpsftP68nLUn4Gik+iJzWtX0GTTqwaO3MmYq8\nUffUnTGMtfHUIJz0IKJYrjKpjik14Cctv8xOO2zMkmfS+eOkAqNpjzvdlqjGnlu+y55bbktUlJVi\n+a7uzcDOWeueALxn0kti514OfBKoAfuAhxBvfbaMHLNWWs1VwO8jGvj1wDFE63478gHejvQ6n0Aa\nvBX4BnAWT/4ZTikNfOq6y/irLz7G6V1SUfD5F/djo4j84CaiSpnC+ZeJIW7UZOBNdUqGvlemYHCH\nXDjHH8c+/gDm/CsTY+CnCKtMwq6rxAiUx6S2if/zZ4vJhZdxWRlBmHjgKaNhcl0ySUKYkwyTegUm\nD2O6t8sw+HolHt5uo7qUcq2MxsY+Dlz6ofVRPdbLjcsFthOHMF2bEyMWT1TgjbK77U4PUvHPZ9fu\n9sbJj5iMB6U4z9gHE70xDt0oT1sHgkQSAnmdnjxghqzhtWt/nkyyr8m4QTuFxDuO64WETvfOJjKE\n9y5TOfIzBt14ycob6LSmnB5A473SIJv63rzG7fV0P2ox1ZFhxCNvuHak66j47yU+T2bmvrOP6X8X\nIKknMysW4WUUY+S7sr5jc4N60tkos+Mf/nU6PzyWfBJvPtj0NDh5O2f3vjK3+FbA6Z+vznW+nUgW\nnpdQdgC3gB+mynOA/cgs9bcCH3fbfQD4CvDZ+c63lnng3hC/G/g08Gakl3m1W3+PW38PiS50ykso\n139xv5tBx7KtZCntOBMThJjOPpk3cnIUdv8UPHE/0Q9uoDY2QuX4MQobNnHwzo+x8y1/AuUJpg8+\nRunMCfGUOvtdCViXgXDkQfHK+06TSoKNmnjgtiH1uH0ec/mEGJt8T3JBZ/KYMI8dOyCeOWCPPwSd\nQ2K8qxNS3tXlQNvatHju0yNSTXD6WJKdEF+QQJgVbdtnlwCmNCjGI8zIkOq4NrbzsPzFb0Kouxzq\netWlCvrJALyRdy5YHEzzxscmKYPY5DYdxNiC3JnUy4msEtc78TIL8t14fTyqujkunYEM83KsbFae\new+VlNGnkOzrvfu03lubcvVj/ACqIDF43rCnA7K+4647+cNn3vgOKs7uwXUs1STYCIkhjUjWNxJ5\nTL4j/91mkjox6TuJWP9G/gv+jizdbiPykcm5eU19+p9xspGPkaRoRXrebFa4CR9E9O3PA69CdPIX\nzLPtgjZwrQz4N90CcBy4Zp7t3uWWpwS2OhmXiN2QhwuueBrBwBbqBx8hs+Us7BMPYy68Bh6+DepV\ngme/kvymcyn0n4WtjrPzoZtgeC90DVB65jWSfXLhC+GhW2Xi33wpmYShPA6jB6Wmd21ajF+hTy7S\neB5J59H67BNncEXzLiRGouc08aAnD0uqoLttt/WyaOTVSTFA+S5MaUNsvE2xD6yVYfJ+NKbPqJgr\nf9hrqSaIZZgZenPsfTVm3qb7GiAzhmRnEsNjvCF1Bs2XQfXHSw+Z9wFH70X6Cn8GVwzK68BO2nHf\ngdQ2qaQ6IDvT265OJKVsjQ8F+awSk2jw3tOE1HeSmoWmXk88YR/IjdMDGyJl+TZ4fbsyKsfzcQP/\n3QQZMZZ+2rsZFRQhHlzjgoeLTSjcerO7cszXiXx3OOLW4SWnK15GYgM/g3jasAQZ2aMjMVvIv734\nDLKBoT8vkzE0JifIjBwms2EbVKcxp18kF+vmc+HHX4Xn/YrUEKlNiiRxzs/A1KeSioI7LhDjvWEH\nM+pVHHtUZtABuSAL3c6LdR5nw02DVhoUw+snEDYh2MkkjzjMSuaHNzb5bslOsQ1MoRdbq0n516iO\ndQWobKMiNUt8fQrA5EX+s7VJkVQaPmcaOa4xkmlSSOq0mLRs4XLPZwTC0vnLtXSGRTTDoEmtE++B\nI/v4EYg+ja9RS/af0QFAPL3c7EEkXgPPFpLOw1o3otN70OWkjSaQqoyNqvyW3vv2QdH4GKlMkHTq\npf8egozUifGedK5T7q5MKINLVuev+9Rjni/yORsDnrMx+S/8zb2NuTecyUOIrPxN4PnIoEaALyIy\n8nsRGfls4HsLHUgNeAt5dMLQlYXeHGzfvZWoVqU2OkKuezApE3pkLwxsg5e+Q2p/VMZkirGoIRkc\nz/hF7Hf+D/RtFsPTNSD1uyeGxVBveRqMPSHrfOGn6qTIIROHpZiVCaF7q5uxBqfrOv3Ua45hVrYv\n9Mj5Rx6Bjo2Yzs3x7brJd4u3DtLW2pRo4Y0aJnSpga44FZmCGO+ojrVWPosvpNSoudopzmutl8UD\nL/ZJDMDkk5TEXCkJvtlaUqPEz76DRWZrSeUb+5lv/IhCiDMcrNPybW3KzZkZYitT0u56EgQ0GXfn\nAYnGbCMJ7HqpYGrYpWFWYOygTJIR1aWcQc9WKPTKHUq+M6m06NrTDrKBknASv8cnEWM9iOjcfwL8\nKvCPSOGJafcaliEjqwFvMbkANuQttYkJChs2ktt6BhQ7k6j/0NmwaVfivRb7XV0RGdloin2SZVKd\nhKP74bQL4dCDMq/liSfEkFWmXEEqI4a91C9eZ6mfeJDG6H7ZJtcl521UZXh9dcLVAkkCbXbiMHRu\nFkNXn5b63VPHkjodTjfH1LBWpATr9Gtbm0o861jzdWmGGIwfwGMCKSXrb+d94avquMQIvMRTr0rQ\n1LhslyCMjbeN6pi4IxLvNZ5b03vITie35RPyuU0owVonIdiG6Lh2akLa67Ju7LEnXFphlAQNOwZE\nfij2y+/SMSifz4Sw+RIJAKthXpcEy69U9bp51j9rnvVLkpHVgLeI6OCPKDegPw9dOSgObSG3YQt2\n/JhMeQbO6z0Gu05ztUhEWzVOHjGFXpEfTn8u3PF56OyDo4/ClvNENunyendWgpfZInRvToyOv2X3\nqXDeW/RTm/kAlJ9goNgnj9PHJXA5cRiyJck2Sde9tpF4216KcfNdmmwHtjImaYLZDnA5zCZTjA09\nxsjdhfdIrZXXQRZbHkuyNUI3G4/Ls7ZuIgVbnUzymF2HASQyBDi9vpIEACtjkqFz5EG5Iyn1QJjD\n+logQRZ6d8hMQH4y5E0XqDF+KtGmv7Ua8FZx9w2UGzJwp6uvJMagewAzuA2OHYTejWKA+7fLwJhi\nn6slEsSRfMHKpL5bz5ec79HDMHJA9i1PyLyLJZdy2unKv9amxUhPHRPpxBtCSGSG8SfE0HcNYTo2\niVfqBkeY7m3YyWHxKHPO+64ed7PJS0aICbPi7WZLcf61nTgknZAfJFSbhCAbe8WAHKvYhy9IZCtj\n4gFHdaluWB2XbXzgszaVZLj4zmL6uHwGE8DkUXn/2KNSPqA8CV2DcPplMLAL070lDlaac9vzIlVa\nT5vabzXgLWP8KAZ4fMpQGC8zmMuLcevZKAa3f7tkGRT7iGtdxNkIPsBWE628OglbL4Xhz8ht/eAO\nGUI/ckDSCX0ZUmNg8piUiQ3zMsGmz7f26W1hHqaOQv8Z+Ip3tjySpIjlOrCjj0rAExLNu7RBjGmj\nCrXJJCjn85jB5Zpn3AhIJ3VkXObL1LBkw5gAxg5go7pkyQShBE7DHHbkYXm/PCqB1sq46MpTx2Hv\nHdDZK9LR5gslUBrmYPtz5I5FUU6Cdr3b0n92q2jUyYdw9kDIwLY+onqNMJuH4Udh63kSrPM1LzIF\nkQYyfsKB1DyOdZfjG0aS693jJhUeOyzbH7xbSsp62aR3u3jeBZdt0aikRlZOJ0PjfZZFZTw1oMan\n8JVku1yHzMZTHEjymKN6UuHQz1RDI6kjXZ1MZBqfF10Zk/X1cmoWGydtjB2UAGx1Cg49IG08/gT0\nbICeTTB0LmbXC+CK3wZaU9JTOfVp17+VGvBWYSO6spDJBtQrVTp7B10Nky6pY3L8AGw6WyoD1iYl\neBk1sJVRV7o1m+Qtu4wU27dNhtNnsjJt2tQxkQr8tlEdpqdk6Hx5NBn04gJzdGxMDLqfPzKqQbUm\nRro6IfJEkIWcm3W+0Ct5xX4IfqMqMoXvFEDe87JH1JB9Jofl9fRIktI3ekgM+RMPQNTAHnkUky/B\nz/wOZudVcP5/WjT3WFFWhTa14GrAW4gBil15Spu3iOH2oyj9LX9WUu0kQ0IGeUjOsMU0JMPDuFlb\nbKMGg2fDvlvFeJcnknzhkQPQuzUZWTj6uKQY5ktJjnTnkBSr8p6zH7QzcQR6d7giTY1kMEw8Gq+c\nGObp49LObElGdTZGpLPIdsg5Jo/KdhNHxFgfeogTt3wNY6Dn6pfBc9+M6RqC5/SpJ620FSeRhbKq\nqAFvFUFILoRapQFRRO2JvWQveQEMni5SQd9W2HCezN5eGoxzp0UmQAZwBPl4FJ4Jc9hCn/OkByS1\nMJMXzTvnRjBOjUCxRzJR/DDoyjgMnOUmJJ50NcFdbZSxAzJMvXwiGZnYqCa1q2uTYpzLo8nozjAL\nww9JHZaxw9KZNBqw7w6m9z9CkM2SP/siuPrXMZf/Fv2vWtmyn4qyGrSrQ6EGvFVYy0DeMjk6TefY\nCTpe9EY4tl+CcpvPleHvtalYOjGFHmx5FEr9kh+drj3hh6qHWeyJQ1IbfPsFkhIXZ6SMieH2s48f\nfwz6TxMDPXU0mRghzCVTW/XukHKy4OSWWjLCszoJx/fD5vPdNGVZGBmWzmL4YSh0Yu+5hanH9jJ1\nYpIN176S0i98EHIdbXsxKMp8tOtfVg14C5moGzafs4XiaWeJdt01KLJG72bxtrOleMZ2O31cBq80\napJf3XAzljeqELl/V5CRqdSsleP5nOaoDuPDsP8OqUx4bB8MnpEMm/cjNH3t6vKYGORsSYKHhW45\n3/SobFvslQ6gUYOjD0nueTaPvftbPH7b7VTLNc54z/8luPw36cp10NXSb1lRVoA2teBqwFuFMWSM\npT41SVSZIvBTdWXdZARdW4hneo8arnY3xKVGw0BGZPocaO+J92yGvd+T14UukTrCAHq3SDna44+J\nvJItiRedLsFaGReD7XOsR/ZJR1CZSEaGlsclwDo1SuORO2hMTnBi3342/szLCd70EU77jf6WfJ2K\nspq0qf1WA94ySj1UI0O92sDWanDiMGw9R/K/p0ck/zvbIduGWai5+SKzyVyNxteZ9h55FEH/mXDw\nJ2LAq9MiacRV9/LyvOCKUE2fIJ4ct2+HnLfULx5654DztrvFez/2GEwcp7b3Hsb37+fQgXF2f+Bm\nMv1nMqR51sopTrvKfnrltYp6lUJo6bvk2dA9INJJtgiH7ofTny350C7n21ZlujHrZ1h3Iwd9cai4\ndkh9GjCSgdK7JfHKa9MSSOzZKh51kIHRR10N60CKPx25T+bLnDwmXvmJgzKM/+hjbtLdKj/65BcI\nDVz0v7/BwIbzxPtXlKcAbWq/1YC3jHwH1chQfuxBCmcVYPA0WT+wI5lZ3E/lBaKFx4NfTBy0jCvu\nYTC+JnfvFjlGdUr09Nq0eNRjT8DATvG8bSQdxomDct5sEQ4/ANWyDAgaPQzD+6kdfYJHv38PZ73+\nF7j08wfcpMSK8tTCzD3fZctRA94qejdL4b7OLhg6U9L7goxUsHOV9ky2JPnd8TB3Nw1YmJrqKj3B\nq6vSZ/MdYpgLXXD/t2DjmeJF0ynyiK/uF4SS9dLhhqwffQymJ6juf5BDdz1Iox5x+vtu4Oy3Xdy2\nt5CKsha0699fDXir2P4MxmtQPvQ4+R/dSJDLEz79GslE6d4alz81QSj1tP0Q83jGmBBws+j4vG0/\ndP3EQfGgM3nIdzgDnYHuDXLuyePioR9+QLJUjh+AE4epHXyEg7ffx+OjdS7/13ul+p6iKG3rwLTn\nfcFTANN/JoGBbG8/2Z27Cc+/Qt6YdtNdRQ1sbQpro2QSBKd/W5ezbX3Z1CAjZVP9zPMDO1w+9l4Z\n8n7iIGDg/htFPpk4KoN6psdl2Prj91N55G4+d8PddPQWuPzfHlPjrShpTJPLGqMeeIswhR46s1Af\nGyU/flwmLp4ckcE2fgJaN3O5RaoJ2npZaoFENWwjwBg/q0wkEw9XxmSy2J6t8JMbZVi+n2vy4E8k\n++TYfgleTooObg/v4/s3fJdN3Rlee+OjqZlhFEXxqAauPInQWGyjDiU31GVqRApCWZnQwNamXODS\nyEw0uS5sVJPsFBsl1QpdlUCT63Lri1CrSnEsEEM+etjVIzlB/fGHCYslDv/wdr56f4U3/vW7CC77\n5VZ9DYrS/qiEoszm4j7L+PEpmByT1L+OAZloOKpLNgnGBTGNzB/p5BQT5tysPG7GHIBGFVsdd380\nC/miaN4Dp8msPbWynGP8OJnObvZ981aeODLNL33lDjXeirIYxjS3PJkPAYeBu1LrrgMOALe75UWp\n994BPAjcB7xwsWapAW8h2WxAsSuPLU/A3h9Iip+b1MC4kZcmzGJrLr/bT/yLlQkTgowYeBuJkc8U\nZaqwbEkCl9NjMgBn3M1KMz5CbWSYh2+8mZGpiIt/8zcxnZta+RUoyrrABGFTyxx8GLh21jqLzDx/\nsVu+4tbvBl7jHq8F3s8iNloNeAs5/S/+heNPjMnoyg074cCPpW62tZJ5ku+Sx6yUjCVy+eEY0cHr\nZRngU5uUgCYWE2Qw2Q6ivXdI6mC+Q3Ty6jTkChy7/2F+eBQu/fzjBC/8kxZ+ekVZRyzfA78ZGJnr\niHOsezkyi30N2Ac8BFy2ULNUA28hwdnXSnpSmJVl41nQs10Mro0kKJnviSc8kNnXJQtFZn0PZLi7\nn+QgzIm3HtUIzrtccrxH9kGjTnnvvUwcGubo0Wlec8N9OopSUZbAKqQR/hbwi8APgN8DTgBbgFtT\n2xwAti50EDXgLSYIAxn5OD4Mm3a5wlSh1Pc2AbZRkWyTbAFTr7oaKPlkNh4bJEFPNx2ayXdjs3nx\nuvNFanvv5oEf7eV7xwJ+5RsPYgo9Lf7UirLOmGcmqG/tL3Pz/vJSj/ZPwJ+75/8T+GvgzfNsaxc6\nkBrwFjM5VWfizlvpfE6P5HBPHoGuLa6uCTJgJ1fATg6LN+4nVbARECa3bUEogU1jZAaf44+7eSWn\nqR4b5oFxw5tefanMCK8oypIw88zIc9WOIlftKMav//LWsWYOdyT1/APAl9zzx4Htqfe2uXXzspoa\neAG4DfgxcA/wl259P3Aj8ADwdaA3tc+SIrCnAjuftpVMZ7ek+nUMytyTDT9RcU6CmW6iYFubTEZh\n+ronjZrkhter2OljMku8CaB3CCrTUJ3myN5hfu5NLyX81Y+37YgyRWlrlq+Bz8Xm1PNXkmSofBF4\nLZADTgfOBr630IFW0wMvAz8FTLnz3AI8D3gZYsCvB/4QeLtb0hHYrcA3gF1AtIptbDmliy7n6H/c\nQGFsWPTsqO4mE65LBcEgi40aUC9jCn3yJ3HlYU2Yi0dsEmRd+VkrnnejDlOjTO17gPFyhPm596p0\noijLZfmTaX8SuAoYBPYDfwpcDTwdkUf2Av/VbXsP8Gn3WAd+nRZLKG4kCTkgRKKxL0M+EMBHgD2I\nAZ8vApsW9U89Nu+i0HUT9G2RgTl+ooVMAbCSaZLrgFzJGetyMglDvZLkgbssFJwOzuhh6OrnG987\nzEtecrEab0U5CU5iJObr5lj3oQW2f5dbmqKZVr0XOL/ZA85x/B8jiew3AXcDm9xr3KNPRN6CRF09\ni0ZgTwmGziXb0yulX8cPi44dhIDFFPpkLszatHjU1kqVwozMVk+Yw+S7MaVBMfiBy2apTUGhkydu\n/CqTdUPmLf+v1Z9SUdY3KyuhrBjNGPB7gf+DaDG/BizFlYuQW4VtwJWIpJImVQ91Tha8fTglyHWQ\n37xDJiIe3CUatk8RdTPAy6Aemxp56Qy5q81tyydEB7cNyRPPd8PhfUSNiJ99pisVqyjKsjEmaGpZ\na5qRUP7ZLecCb0QE91vcupuaPM8o8GXgUsTrHgIOIWK+j8g2HYG97p3Xx8+vvuJyrr7y8iab0X4E\np11OtOuHcOQRmDiE2XCevBHmUrPpVAELQRaT705qgod57PQIJtshAU5rwboKhX2bGBursvV//G1L\n/liK0gr2fOvb7Ln52yt/4DYN/jergYeIAT8PGAbuAH4X8chfM88+g4gQfwIoAi8A/gyJtL4BeI97\n/ILb/ovAJxDJZisLRGCv++O3NdnsdcTwfghCySwJcxgbYSvjmCCLtXUpahXVXfqgEUmlIV64rYyJ\nhNKoyHvTI3B4H6EBM3RRqz+ZoqwZV18506H7s3f91Yocd740wlbTjAF/H/CzwH8A7yQxqu8B7l9g\nv81IkDJwy0eBf0eKt3waSVzfB7zabb/kCOwpw/aL4f7vQr5bvOnI1TexUvPbBNlkYuIgI4WrooZI\nKLYhOrnPYLEWalNElTLd/UWtdaIoK0Gb3sU2Y8DvBP47MDnHe89aYL+7gEvmWH8cuGaefZYUgT1l\n6NjI9KMPUqyXsY2qSCb5HCbfJSmBtoFtVJNZ6DGpuSmNGO/6tFQrzHdjGxWiaplMTrVvRVkR1rGE\ncidwzqx1o8CjiDyinCRm4GxJU6qMSwaJqwEeSyaZgisf20jqgDsPXQpYhdhMQTz2mmRu1icn6Nq6\necHzKorSHOt5Qod/RIKPd7rXFyDpgD3AfwO+tjpNe+pgTEDhgueITJLJA0gNlEyBOCPFTdoAiIce\nNdw8mDIvpjEhtiG1xJkeYeLwMIN//pW5T6goytJoUwmlmVYdRFIBL3XL04FHkKDk9QvspyyFZ7wK\nHvw2JsxjwqwYbxPIrZv3xI1x2Sc5Ga1pIzesPgtY2acyBtki3du3YfrPbPWnUpRTAmNMU8ta04wB\nPwfxuD33IBkpD/NUCTKuAWbDebD1fGx13AUkG7M2cLPRW5sENKMGhFlsoxLXTiFThGOPyPB7RVFW\nhsA0t6wxzUgodyPlDz+F3M+/GjHieWTYu7ICmNIAtjwupWQDKVYVyyahy0LJFp3MkkvklrorN4tN\nZuPp30n+0mtoyTTZinIK0q5jKZpp1RsQb/u3gbci8skbEOP9/NVr2lOQM56HKfaJt22TGl62PCYG\n21qZtKEyDoGTVzI5KT0b1aE8CrVJqE7IrPOKoqwMbTqUfjEPPAPcgAyB/19zvD++4i16KhPVsccf\nxnRvc3q3ZJuYQo8Mk3cGXF7bJKWw2AdRA9OxEVsZlf2e8aq2DbwoynqjXbNQFmtVHaln0rvIdspK\nUJ2A0oCTTdwISxvJXJhRQ9ZXxmWG+RlZKfVEM48aMsHxwC6t/a0oK4UJmlvWmGY08ElkUM6NJIN5\nLPCW1WrUU5ae7WKgsyVsdTIZyIPLPoka8ugNdmDlPffHsbWppBxtQWfeUZQVo02doWYM+Ofc4jNO\nDJp9siqY0gbssQexxX5MriRSSGNCHo0B4+qggAQ0a1PExr3uRmpisc6DVxRlZWjXu9lmDPi/ACXg\nNGSqM2WVMIVubL5L/izWuqnVMon+jZHXGfe+v2WzVgx8kJUAZ9/pbRs1V5R1yTrVwEFm0Lkd+Kp7\nfTFSOVBZDcKcGOl6hbhcepCRzBSspBRmnEduo3jyB0CmXQtzaPqgoqwwy89C+RBSQvuu1Lq/QuZZ\nuANRN9JzLCxpXuBmDPh1SNGqEff6duCMJvZTloHp3ekmMS466SSEepnYSNtIhs+DDNoxYZydAmBr\nk1qBUFFWmuUHMT8MXDtr3deRWc4uQiZ3f4dbn54X+Frg/Sxio5sx4DWeXLTqlJ5ouJWYjg3yJKq7\nGWceSt0AAA7uSURBVOjdQB6TkaVRn1HISmQVV3UwzGHyPap/K8pKE4TNLU/mZhLn13MjiQ29DZm8\nBuafF3j+ZjXR9LuBn0f08rOBvwe+08R+ynLJdYqRzpacdGKSgTveGzfGed8NVwsFJ7U4g68oysqx\negN53oSMtYFlzAvcjAH/LcTdryC9wxgyKlNZLbyObQJolJ2MUpEePlsU7zyencc6o46TXbJSelZR\nlJVjdfLA/xioIjORzceCGX/N5oH/kVuUNcDkOuRXq02Lzh013Ew8NTHmuY5ESjGhpBHahhhyo8Zb\nUVacebzrPQ8c55sPHl/OEd8IvBj46dS6pucF9jRztZ8D/D6wM7W9ReugrCom14EtnwAbuImKcyk9\n3N+uuZno/XMkFVFRlBVmHu/66nMGufqcwfj1n9/wcDNHuxb4A+AqoJxa3/S8wJ5mDPi/ItUIPwD4\nGqU6kGctCLIyMjPXKfKJjaA6Ll45OMMdymTGJhDNXFGUlWf5A3k+iRjqQWA/8KdI1kkOCWYCfBeZ\nA3jJ8wI306ofIhM5tAvWTg63ug1rhq1OJnVQ/AAfV0Y2/m1NINv4QUCKogBxVtfJXhS28f+9pKkN\nw1/78kqcr2maUd2/BPwGMst8f2pR1gCT63Dedz0ZtGMb8hi4DJWoDrmSGm9FWS2Wn0a4qjQjobwR\ncfV+f9b601e8NcqcmFwH4LxxXPqgaUggU71uRVl92rQ0RTMGfOdqN0JpDm/IAc31VpS1pE2dpIW6\nlbelnr9q1nvvWoW2KIqitCdtWg98oTO+LvV8dg74i1ahLYqiKO3JOp1STVEURVnHGriiKMpTm3Wo\ngV+ITFo8DlyQeu5fN8N24CakINZPSKZh60eS2B9ASium59xcUj1cRVGUVadN0wgXMuAh0OWWTOq5\nf90MNeB3kGJYz0byyc8D3o4Y8F3Av7vXsIx6uIqiKKvOOgxirgSHgB+75xPILBRbkVl+PuLWfwR4\nhXu+5Hq4iqIoq06bBjHXssvYiUzHdhuwCZlmCPfop5BZcj1cRVGUVadNPfC1CmJ2Ap8F3opo6Gn8\n1DLzoYWzFEVpLW0axFwLA55FjPdHgS+4dYeBIURi2Qwcceubqod73Tuvj59ffcXlXH3l5SveaEVR\n1h97vvVt9tz87ZU/cJumEa52t2IQjfsYEsz0XO/WvQcJYPa6x91IPdzLEOnkG8BZzPTCn1LVCBVF\nWT4rVo3wU29uasPwtR9cifM1zWp74JcDrwfuRGazB0kTfDdS9/bNSLDy1e69JdfDVRRFWXXCtU8R\nbIbVNuC3MH+g9Jp51r8LrbWiKEo78RTWwBVFUdY3baqBt2erFEVR2omTywN/K3AXMhr9rW7dQqPR\nm0YNuKIoymIsPw/8acAvA88ELgJeCpzJ/KPRl4QacEVRlMVYvgE/Fxm8WEYmhf8m8J+ZfzT6klAD\nriiKshjLL2b1E+AKRDIpAS9GxrfMNxp9SWgQU1EUZTGWH8S8Dxnv8nVgEqkN1Zi1zWKj0edFDbii\nKMpizGPA99z5GN+8a/9ie3/ILQDvRGo8zTcafWnNWs5OLUZHYiqK0hQrNhLzy29bfCsgfMn1c51v\nI2KgTwO+hpTW/mPmHo2+JNQDVxRFWYyTywP/DDCAlMn+dWCU+UejLwk14IqiKItxcgb8yjnWHWf+\n0ehNowZcURRlMdp0JKYacEVRlMVowXyXzaAGXFEUZTHUA1cURVmnBGrAFUVR1idaTlZRFGWdohKK\noijKOkUNuKIoyjpFJRRFUZR1StCeprI9W6UoitJOqISiKIqyTlEJRVEUZZ2iHriiKMo6RQ24oijK\nOkUlFEVRlHWKeuCKoijrFE0jVBRFWacE7SmhrPZ9wYeQyTvvSq3rB24EHkBmau5NvfcO4EFkJucX\nrnLbFEVRmsMEzS1rzGqf8cPAtbPWvR0x4LuAfyeZyHM38Br3eC3w/jVon6IoyuKcnAHvRebFvBe4\nB3gWCzuyTbPaBvJmYGTWupcBH3HPPwK8wj1/OfBJZOLPfcBDwGWr3D5FUZTFOTkD/rfADcB5wIWI\nwjCfI7skWuHhbkJkFdzjJvd8C3Agtd0BYOsatktRFGVujGlueTI9wBWInAxQR2aln8+RXRKtliis\nWxZ6X1EUpcWYJpcncTowjMjJPwL+Gehgfkd2SbQiC+UwMAQcAjYDR9z6x4Htqe22uXVP4rp3Xh8/\nv/qKy7n6ystXpaGKoqwv9nzr2+y5+dsrf+B5JjXe8/17+Ob371lozwxwCfCbwPeBv+HJcslijuy8\nrEVuzE7gS8AF7vX1wDHgPcgH6XWPu4FPILr3VuAbwFk8+YNZOzm86o1WFGX9Yzo2wMnbOdu46/82\ntWF4wWtmn28I+C7iiQM8D8m2OwP4KRJH9ibg3KU2bLUllE8C3wHOAfYDvwS8G3gBEn19vnsNEp39\ntHv8CvDrqISiKEo7sHwN/BBi+3a519cAdyNO7RvcujcAX1hWs5azU4tRD1xRlKZYMQ/87n9tasPw\n/FfNdb6LgA8AOeBhxJENEYf1NCTr7tXAiaU2TEdiKoqiLMpJiRV3AM+cY/01J3NQUAOuKIqyOFqN\nUFEUZZ1i5s5CaTVqwBVFURZDPXBFUZR1itYDVxRFWa+oB64oirI+UQlFURRlnaISiqIoyvrEqAeu\nKIqyTtE0QkVRlHWKeuCKoijrFdXAFUVR1ifqgSuKoqxT1IAriqKsV1RCURRFWZ+oB64oirJO0TRC\nRVGUdUqbeuDtKewoiqK0EyZobnkyBeA24MfIfL9/6db3AzcicwN/HZncfcmoAVcURVmM5U9qXEZm\nn386cKF7/jzg7YgB3wX8u3u9ZNSAK4qiLIppcpmTKfeYQyYzHgFeBnzErf8I8IrltEoNuKIoymIs\nX0IBsbM/Bg4DNwF3A5vca9zjpuU0S4OYiqIoi3JSQcwIkVB6gK8hMkoa65YlowZcURRlMebxrvd8\n5/vs+c4Pmj3KKPBl4FLE6x4CDgGbgSPLatZydmox1k4Ot7oNiqKsA0zHBjh5O2ejQ3c2tWEwdOHs\n8w0CdeAEUEQ88D8DfgY4BrwHCWD2soxApnrgiqIoq8dmJEgZuOWjSNbJ7cCngTcD+4BXL+fgasAV\nRVEWY/kDee4CLplj/XHgmmW3x9GOWSjXAvcBDwJ/2OK2KIqicJJphKtGuxnwEPgHxIjvBl4HnNfS\nFi3Cnm99u9VNmIG2Z2G0PQuj7ZmH5Q/kWVXazYBfBjyEaEI14FPAy1vZoMXYc3Ob/MEc2p6F0fYs\njLZnPtQDb4atwP7U6wNunaIoSus4uYE8q0a7BTGXlcyuKIqyqrRpNcJ2a9WzgesQDRzgHcgopvek\ntnkIOHNtm6UoyjrlDmQU5Mlgo6MPNLVhMLgL1tCutpsH/gPgbGAncBB4DRLITHPWGrdJUZSnOm3q\ngbebAa8Dv4mMVgqBDwL3trRFiqIobSdWCO3ZKkVRlPbBRsceamrDYOAsWEO72m5ZKIvRikE+H0IK\nz9yVWrfQbBrvcO27D3jhKrRnO0lJyp8Ab2lhm5Yz28hqfz8gd2+3A19qk/bsA+50bfpeG7SpF/gM\ncnd7D/CsFrbnHOR78cso8p9u9W82A2OCphZlfkIkgLkTyCJGYy0G+VwBXMxMA3498Db3/A+Bd7vn\nu127sq6dD7HyneQQSVCmE7gf+R5a1aaSe8wAtyKzjbTy+wH4XeDjwBfd61a3Zy9ikNK0sk0fAd7k\nnmeQMqet/o5wx30CcVLaoT0ea0f2NrWgmXTz8hzgq6nXb2eZ0xAtg53MNOD3kRRgH3KvQTyD9J3B\nV5HMmtXkC0hNhVa3qQR8Hzi/xW3ZBnwDqbnsPfBWfzd7gYFZ61rVph7gkTnWt/o7AvGmb26j9nis\nHdnX1MIaG/D15PO30yCf+WbT2IK0y7PabdyJ3B3c1sI2LWW2kbX4ft4H/AGSfupp9e9lkU7lB8Cv\ntLhNpwPDwIeBHwH/DHS0sD1pXgt80j1vh/Yk6FD6k6Zdb00W63VXq92dwGeBtwLjLWyTn21kG3Al\nS59tZCXb8lKkMP7tzB9IasXvdTnS0b4I+A1ElmtVmzJIdbz3u8dJnnwn24rvKAf8LPCv85yvFddY\nCh1Kf7I8jmhjnu3M7InXEj+bBsycTWN2G7e5dStNFjHeH0UklHZo01yzjax1W56LTBa7F/Hkno98\nR63+bp5wj8PA55GaP61q0wG3fN+9/gxiyA+1qD2eFwE/RL4jaP1vNhP1wE+a9CCfHDLI54sL7bCK\n/P/t3T9rFFEUhvEnGgUJqfxTWJlCG1ELbQRBC4WIYKE2FhL8EDZptLGzELS0FOwsYqGFiPgFjIii\noG20sNJCCOhYvLPcmASzWbO5O+H5wbKzYWEPw+bunXvPOTMHzLTHM5RBdI5cBu4kl6oHKVkHG2WM\n5Me/B+5WjmkPJTtgF3COzH5rnZ9Z8s891X7OC+BaxXggewOT7fEEWed9WzGmr2Qp8lD7+ixZ9npS\nKZ6eq5Tlk97n1oxnmdGcgXfNeZJ18YlsZmyGR6QqdJF88a+TjILnrJ7iNNvG94HcNmmjnSLLFvOU\n1KvpSjEdIeuo8yRN7kb795rnp+c05Qe+ZjxT5PzMk7TP3ve2ZkzHyAz8DfCYbGzWjGcC+Eb5oaNy\nPMs1zfeFvh5s8lKvPxmS9G9N82OhrzeOTe6HlePqNLla3g484O/eTv+lS0soklTJwEsoQ71JjQO4\nJK1l8E3Mod6kxgFcktY08Ax8qPUro9aNUJJGz+ApgkPd1HQAl6Q1jE3s6/ety4vqRql+RZK0DuPA\nZ0r9ymY14ZOq2k3Jc/9CZi2vyQznfsW4pPWqUb8ijYybpEWspCXMQlFX9HaRzlDaxN4iva1fkTSt\nS8AdUhX6lLLHcxx4SdoxPKP02JA6zQFcXTdFOiBeBB6Su7gcBX4CF0jjr3vAZeAEaaN6u0qk0gYz\nC0Vd1pCZ9i/SZ2QbuSE2pGHUAdK06TDpqwGpjOuvLloacQ7g6rrF9vk3qXRjyetxsvTyjrSalbYU\nl1DUZf1UV3wE9lJuu7WD9KSQOs8BXF3RLHle7RhWVr01ZFZ+hXSA67XgPTm8MCVJkiRJkiRJkiRJ\nkiRJkiRJkiRJW8oficN/jVry+kIAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "transfer.plot(response='2d')" ] }, { "cell_type": "code", "execution_count": 48, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEZCAYAAACXRVJOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcU/W9//HXzLAN67DIjgKj7KKgBUTQoGLBFXGrC2Kr\n9v5aH9Vae4va3jpqr4pVq9dWbd2KtLV1qV6xoqKXWMRdQdlRFtlBZZedye+Pz4lzCJmZzExOvifJ\n+/l45DEnJ8nJZwJzPvl8twMiIiIiIiIiIiIiIiIiIiIiIiIiIiIiaTcXOMF1EGkUAVYGcNyuQDlQ\nGMCxRQJRz3UAEjrbgZi33QTYBez37v8Q6OciqBy0HGiLfbbfANOAq4GtDmMSEanUMuAk10F4gvry\nEsFNheD/bNsBs4G7AohDJGUqZ6WmllNxIisDngEmY99sPwWOAG4E1gNfACN9r20BPAasAVYBt1H1\n/8Fy4MfAZ8Aib98Z2MlzEzATONL3/AnecbcCC31xNgTuA1Z7t98BDZK83wTv9/G737tVF38hcDfw\nJbAEOL2K3yvReuA1oK9v3xDgbez3nA2c6Hvscu89tgJLgYt9+2cCDwCbgQUcmNA7Ai8CX2Of6ZW+\nx8qAp4FJ3nHnAsf4Hq/ssy0AbgA+B74C/gG0TPH3FpEskqxC8O8rA3ZiJ/0i7GSyHEsIRdgJZ6nv\ntc8DDwHFwCHAe1gzVGXKgVeBEuykPgA7eX4HOxFd5sVTH+gJrADae689FOjubd+KnVzbeLeZ3j44\nsEI4DGu+aerdL8JO/oNSiP//YSfgTtgJcTrWHFRVhXCyt90ZS6a/9u53wk6uo7z7p3j3W2PNeFuw\nxAtWXfTxti8H9gLXerFfgCWGEu/xfwO/x5LhUcAGYIT3WBn2bzkK+2xvB97xHqvqs70W+2w7Yv8O\nDwN/q+R3FpEslkpCeNX32JnANuyEAtAMO6k3x05cu4BGvudfBPxfFe9fjp2w4x6i4kQetxDr5C7F\nksXJ2InJ73MqTq4Ap3q/BxzcZDQDGOdtj/ReSwrx/x8HJreRVN1ktBz7rLZ6z3ve99wJwJMJz38F\nS4CNsaphLJaY/C7HKiC/94BLgS7APiyhxN0OPOFtl2FVSlwfYIe3fTiVf7bzOfD/SAdgD2p9yEr6\nR5O62uDb3ol9k4357oN94z4MO5msxU5om7Bvk4d4z5mHnSC3Acf7juk/WR8GXO97/Sbs23UHrAnl\np9iJbT3wlLcf7NvrF77jrPD2JfM37EQP1hTzV997VxV/h4RYV1Ry/LgYcDaWLCPYSfVY33udn/B7\nHo99Q98BXIhVJGuAl7Bv8HGJCeELL7YOwEasAvLH2Ml3f71veweW/AqxpFjZZ9sVS2bxOOdjiadd\nNb+/hJASgmTKSmA31uzR0ru1oKIPoC9WUTTDmnTiYr7tFcB/+17fEks2//AefwoYjp1QY8BEb/8a\n7MQVd6i3L5lnsRN0J2AMFc0f1cW/1juu/z1S9W+s3T8e7wqsX8b/ezajotP5NazKaY9VSI/4juU/\nwYN9Fmu8WysqmsPiMa5KMcbKPtsVWPXlj7Ux9nlIllFCkExZi53I7sVOboVYM09N5jQ8gn0zHoQ1\nSzXBOm+bAj2wb9kNsRO3f7jsU8CvqOhD+DV2wk3mSyAK/Bnr/4h3ZlcX/9PANVT0IdxQg98LrNN7\nEDAY+AvW/HYq1hfQiIok1RarLJpg/QXf+H5PvMevwaqZ84FewMvYif9t4A7sM+oP/MB7r+pU9dk+\njDU9xRPgIcBZqf/aEiZKCFIXMQ78Bk819y/DOjTnY80Xz1DRUVnZ8f0+Aq7COkY3YiNlLvMea4id\n7L7ETt5tsM5tgN8AH2Idt59627+p4n3+hrWXJ3aOVhX/I1h/yife8Z9LctyqfIV1ysdH85wN3IQ1\nya3AmsoKsL/Z67Cmoa+xb+0/8h3nPazD+UtsFNS5WFMOWFNYV6xa+CeWGON9IFX9W1b12d6PjVx6\nDesPeYeKTniRgxQBs4Ap3v1W2CScxdh/opJKXiciNXM51ikuUiuZqBCuxb5Rxb9t3IAlhB7AG9S8\ntBYRkSzUGXgdG+scrxAWUjECId4pJiJ1Nx7roBYJpWewyUQnUpEQNvkeL0i4LyIijgTZZHQG1iE2\ni4qJSomSdWSJiIgDQa52OhQbfnYaNmyuOTbUbz3WVLQOm9yyIdmLS0tLY0uWLAkwPBGRnLQEm11e\nY0FWCDdh0+W7Ad/DhreNw4aojfeeMx54IdmLlyxZQiwWC/3t5ptvdh5DLsSoOBVn2G/ZEic2P6ZW\nMjkPId40dCe2zstibLLLnRmMQUREKpGpC+S86d3AJvSckqH3FRGRFGmmch1FIhHXIVQrG2IExZlu\nijO9siXOuqhs9E8YxLz2MBERSVFBQQHU8tyuCkFERAAlBBER8SghiIgIoIQgIiIeJQQREQGUEERE\nxKOEICIigBKCiIh4lBBERARQQhAREY8SgoiIAEoIeWXqVLjvPtdRiEhYaXG7PDJwIMyaBfpYRXKX\nFreTlDRq5DoCEQkzJYQ80rCh6whEJMyUEPJIvEJQk5GIJKOEkEf27LGfu3a5jUNEwinohNAIeA+Y\nDcwH7vD2lwGrgFnebVTAcQiwcaP93LbNbRwiEk71Aj7+LmAEsMN7r7eAYUAMuNe7SYbs2GE/t2+H\ntm3dxiIi4ZOJJiPvNEQDoAjY5N0P85DXnLR7NzRvrgpBRJLLREIoxJqM1gPTgXne/p8AnwCPASUZ\niCPv7dkDrVsrIYhIcplICOXA0UBn4AQgAjwEdPP2rwXuyUAceW/3biUEEalc0H0IfluAfwHHAlHf\n/keBKcleUFZW9u12JBIhEokEFlw+UEIQyT3RaJRoNJqWYwXdjt8G2AdsBoqBV4FbsGajdd5zrgO+\nA1yc8FotXZFmDRrABRfAiBFwxRWuoxGRINRl6YqgK4QOwCSsaaoQmAy8ATyJNRfFgGXAfwQcR94r\nL4e9e6FVKxtlJCKSKOiEMAcYmGT/ZQG/ryTYuxfq14emTeGbb1xHIyJhpJnKeWL3blvLqEkTJQQR\nSU4JIU/s2VORENRkJCLJKCHkid27rVNZFYKIVEYJIU+oyUhEqqOEkCf27FGFICJVU0LIE6oQRKQ6\nSgh5wt+prIQgIskoIeQJdSqLSHWUEEIkyJU61GQkItVRQgiJdeuge3coKLDtdNuxw5KBEoKIVEYJ\nISR++EO48EIYOBAWLUr/8XfsgOJiJQQRqZwSQgisXg0zZ8Ktt0KfPrBsWfrfY8cOaNzYbjt2BNs8\nJSLZSQkhBN56C4YPt07fbt1g+fL0v0c8IRQVWV/Czp3pfw8RyW5KCCHw5ZfQubNtd+0abIUAajYS\nkeSUEEIg3r4PwVcIoAXuRCQ5JYQQ2Lmz4mQdVIXgfw9VCCKSjBJCCPgrhE6dYO3a9Hf6qslIRKqj\nhBAC/m/vDRpYp++2bel9D3/S0VXTRCQZJYQQ8J+swa57vGlT+t8jnnRKStJ/fBHJfkEmhEbAe8Bs\nYD5wh7e/FTANWAy8BpQEGENW8FcIYAlh48b0vsfmzdC8uW23bp3+44tI9gsyIewCRgBHA/297WHA\nDVhC6AG84d3Pa4kVQsuW6f8Gv2EDtGtn261awddfp/f4IpL9gm4y2uH9bAAUAZuAs4BJ3v5JwJiA\nYwi9TFQI69dXJITWrZUQRORgQSeEQqzJaD0wHZgHtPPu4/1sF3AMoRd0hRCLWYXQtq3dV5ORiCRT\nL+Djl2NNRi2AV7FmI7+Yd0uqrKzs2+1IJEIkEkl7gGEQdIWwZYuNXIonnTZtbHa0iGS/aDRKNBpN\ny7GCTghxW4B/AcdgVUF7YB3QAdhQ2Yv8CSGXJRtllM6E4G8uAujYEdasSd/xRcSdxC/Lt9xyS62P\nFWSTURsqRhAVAyOBWcCLwHhv/3jghQBjyAqJFUK6m4z8zUVg6yatXp2+44tIbgiyQuiAdRoXerfJ\n2KiiWcDTwBXAcuCCAGPICpmuENq1s07lvXuhfv30vY+IZLcgE8IcYGCS/RuBUwJ836yT6QqhqMia\njVasgNLS9L2PiGQ3zVQOgaArhLVrD6wQAPr2hTlz0vceIpL9lBAc27vXhoX6m27SnRAWL4YePQ7c\nd+SRMHdu+t5DRLKfEoJjO3dadVBQULEv3U1GCxZA794H7jv0UFi1Kn3vISLZTwnBscT+A7A1h3bs\nsOqhrvbvh88+g549D9zfrp11NouIxCkhOJbYfwBWLZSU2IJ0dbVsmXUoN2ly4H4lBBFJpITgWLzJ\nKFG6+hGSNReBJYkNlU4JFJF8pITgmP86BX7pWoCusoSgCkFEEikhOLZrFzRqdPD+Dh1suGhdLVyY\nPCE0b259FDt2HPyYiOQnJQTH9u2DekmmB6ZrvaFPP7U5B4kKClQliMiBlBAc278/uISwdatVCMcc\nk/zxdu3UjyAiFZQQHNu3z5aSSJSOhDB9OgwZkrxJCqxjWRWCiMQpIThWVYVQ1xVJp06F0aMrf7xd\nO1i3rm7vISK5QwnBsaAqhFgMXnkFRo2q/Dnt26tCEJEKSgiOBdWHsHAhlJdDnz6VP0cJQUT8lBAc\nq6xCaNnShqR+803tjvvGG3DqqQeukZRITUYi4qeE4FhlFUJBgVUJtZ2LsGgR9OtX9XNUIYiInxKC\nY5VVCFC3ZqMlS6q/+I0qBBHxU0JwrLIKAeqWEJYurT4hqEIQET8lBMeqqxAuvRTmzavZMffvhy++\ngG7dqn6elq8QEb+gE0IXYDowD5gLXOPtLwNWAbO8WxWDI3NbVRXCjTfC4MEwc2bNjrl6ta2WmmwV\nVT8tXyEifpWcitJmL3AdMBtoCnwETANiwL3eLa9VVSG0a2cjhVaurNkxU+k/iItPgKuumhCR3Bd0\nhbAOSwYA24EFQCfvfhUDIvNHVRUC2KUuV6yo2TFrkhC6dbOL6IiIZLIPoSswAHjXu/8T4BPgMaAk\ng3GESlUVAmQmISxdWrPji0huCrrJKK4p8CxwLVYpPATc6j12G3APcEXii8rKyr7djkQiRCKRgMPM\nvOoqhC5datdkdM45qT23tNQmsYlIdopGo0Sj0bQcKxPNNvWBl4CpwH1JHu8KTAGOTNgfi8ViwUYW\nAnfeCZs2wcSJyR/fudNmLe/YYduJ10ZOtGoVjBkDDz4IgwZV//4ffwyXX27XTRCR7FdgyxPU6twe\ndJNRAdYkNJ8Dk0EH3/Y5wJyA4wityi6QE1dcbMNDH3gAmjat+ljbt1tF8dFHqTcZ9e0Ln30Gu3en\nHrOI5Kagm4yOBy4FPsWGlwLcBFwEHI2NNloG/EfAcYTW/v1V9yEAHH00/PSnth2LVb4+kX94aqtW\nqb1/w4bQqZPNW+jRI7XXiEhuCjohvEXyKmRqwO+bNfbtg/r1q37OmDEwbZptb90KLVokf97HH9v1\nD37+86oXtUvUvbt1LCshiOQ3zVR2rLpOZYDx4+Hcc+GII6yPIJnHH4eHH4bzz4eTTqpZDN27W0e0\niOQ3JQTHqht2CtaR/OyzNgS1squoXXGFDU9NpSM5UbxCEJH8poTgWCoVQlynTskrhG3bLGmsXWud\nxDVVWqqEICJKCM6lUiHEde5sFcLChdbJvGmT7Z83D3r1stVLa0MVgohA5iamSSVqUiF07gxXXw33\n3WejjXr0gEsugT/+sfqL4VQl3odQ1QgmEcl9SgiO1aRCOPVUm4F8111WFdx7r406+vOfYWodxm21\naGHDT7/8Etq2rf1xRCS7KSE4VpMKobQUnnvOtjt3ttFHmzdbh/OoOi4gXlpqVYISgkj+Uh+CYzWp\nEPwaNoQRI+CTT+CUU+oex4AB8NZbdT+OiGQvJQTHalIhJLr/fvj97yufqFYT55wDU6bU/Tgikr2U\nEByrbYUAtm7R1VenJ46hQ22m85496TmeiGQfJQTHUlnLKBOaN7d+hNmzq3+uiOQmJQTHysvDkRAA\njjsO3n7bdRQi4ooSgmP790NhSP4Vhg6Fd95xHYWIuBKSU1H+ClOFMHSoKgSRfKaE4FiYKoTSUti1\nq+aX7BSR3JDKqahPkn2RNMeRt8JUIRQUqNlIJJ+lkhCeBiZgl8NsDDwA3BlkUPkkTBUCqNlIJJ+l\ncioaDHQB3gHeB9YCQ4MMKp+EZdhp3HHHqUIQyVepJIR9wE6gGGgELAXKgwwqn4SpyQigf39bOK9c\n/8IieSeVhPA+sAs4FhgOXAw8k+LxuwDTgXnAXOAab38rYBqwGHgNKEk95NwStiajkhJbCmPFCteR\niEimpXIquhL4L2Av1lx0FpDqqjd7geuAvsAQ4GqgN3ADlhB6AG949/NS2CoEsKuuzZ/vOgoRybRU\nEsJ64NCE25spHn8dEF8MYTuwAOiEJZVJ3v5JwJgUj5dzwlYhAPTpY81GIpJfUlln82Ug5m03AroB\ni7Bv/TXRFRgAvAe0wxIN3s92NTxWzghbpzJYhaCRRiL5J5WEkHhxxoFY009NNAWeA64FtiU8FqMi\n4RygrKzs2+1IJEIkEqnh24ZfeXn4KoS+feGRR1xHISKpiEajRKPRtByrtlfQncvBiaIy9YGXgKnA\nfd6+hdjktnVAB6zjuVfC62KxWNI8kVOOOcauiXzssa4jqbB5sy2tvXWrrrEskm0K7I+2Vn+5qVQI\n1/u2C7EKYXWKxy8AHgPmU5EMAF4ExgMTvZ8vpHi8nBPGTuWSElsOe8UKOOww19GISKak0ljRDGvy\naQo0wL7tn53i8Y8HLgVGALO82yhspvNIbNjpSeTxzOcwdiqDdSxrpJFIfkmlQiirw/HfovKkk4Yr\nAWe/MFYIAP36waefwujRriMRkUypKiFUNdcghg0dlToKa4Vw4ol2veYJE1xHIiKZUlVCuKeKx3K/\ntzdDwjjsFOCUU2DcOOtYbt7cdTQikglVJYRlwBeZCiRfhXHYKUDTpnD88TBtGpx7rutoRCQTqjoV\n+Uf+PBd0IPkqrBUCwOmnw7/+5ToKEcmUVL+bdg80ijwW1k5lqEgIWvlUJD+EsLEiv4S1Uxmge3do\n3Ro+/NB1JCKSCVWdivpjy0xsA470bW8DtgYfWn4Ic4UAcMYZajYSyRdVJYQibFJaM6zzuZnvpnEn\naRLmCgEsIbz0kusoRCQTQnwqyg9h7lQGu8byypXw+eeuIxGRoCkhOBb2JqN69eDii2HyZNeRiEjQ\nlBAcC3uTEcB558Hzz7uOQkSCFvJTUe4Le4UAcNxxsGEDLFjgOhIRCZISgmPZUCEUFdmaRuPHw+7d\nrqMRkaCE+fIneXGBnIYNYcsWaNTIdSRVi8Ws6ai0FO66y3U0IlKZulwgRwnBsfr14ZtvoEED15FU\nb80a6NkTNm2yzmYRCZ+6JISQN1bkvrAPO/Xr2BE6dICFC11HIiJBUEJwKBazW9j7EPyOOQY++sh1\nFCIShCw6FeWe8nK7iH02Xcj+mGPg449dRyEiQVBCcCgbhpwmUoUgkruCTgiPA+uBOb59ZcAqYJZ3\nGxVwDKGVDUNOEw0YAJ98YrGLSG4J+nT0BAef8GPAvcAA7/ZKwDGEVjZ1KMeVlED79lYlaE6CSG4J\nOiHMADYl2Z9FrebBCevlM6szciQMHgzXXus6EhFJJ1eno58AnwCPASWOYnAuGysEgGuugf79bT6C\niOQOF9OLHgJu9bZvA+4Brkj2xLKysm+3I5EIkUgk4NAyKxs7lQF69YLHHoMrr3QdiYhEo1Gi0Wha\njpWJppuuwBTsqms1eSznZyp/+SX07g1ffeU6kprbvh3atoVt27IzqYnkqmybqdzBt30OB45AyivZ\nWiEANG0KhxwCy5e7jkRE0iXohPAU8DbQE1gJ/ACYCHyK9SGcCFwXcAyhlY3DTv0GDoRXX3UdhYik\nS5hH++R8k9HKlXatgVWrXEdSO7NmwahR8N570LWr62hEBLKvyUg82TrsNG7AALj5Zhg+XHMSRHJB\nFp+Osl+2Djv1+/GP4fDD4eWXXUciInWlhOBQNncq+116KfzhD7Bnj+tIRKQulBAcyvZO5biLLoLi\nYrjqKteRiEhd6LpXDuVKhdC4MUyaBN27w+bNtt6RiGSfHPh+mr1ypUIAaNUKTj4Znn3WdSQiUls5\ncjrKTrnQqex33nnw4ouuoxCR2lJCcChXmozivvtdiEZhxgz73UQkuyghOJRLTUZgzUZdusAJJ8Ab\nb7iORkRqKodOR9kn1yoEsOrgl7+EqVNdRyIiNaWE4FCuVQhgVcKYMUoIItkox05H2SXXOpXjBg6E\njRth2TLXkYhITSghOJTtaxlVprAQRo/WiCORbJODp6PskasVAthyFpMmuY5CRGpCCcGhXOxUjjvp\nJGs2+uAD15GISKqUEBzKxU7luMJCuOkmazr6619dRyMiqcjR01F2yOUKAeCHP4RnnrFhqPv2uY5G\nRKqjhOBQLlcIcSNG2GJ3LVrA2rWuoxGRquT46SjccrlT2W/GDBg7FiZPdh2JiFQl6ITwOLAemOPb\n1wqYBiwGXgPydrHkXB12mqhZM7jiCpgwAW67zXU0IlKZoE9HTwCjEvbdgCWEHsAb3v28lC8VAtj6\nRtdfD7feCitXuo5GRJIJOiHMADYl7DsLiI9QnwSMCTiG0Mr1TmW/wkK4+24YPx7+8Q/X0YhIMi4a\nLNphzUh4P9s5iCEU8qFTOdGVV8LEifDRR64jEZFEri+hGfNuSZWVlX27HYlEiEQiwUeUQflUIcQN\nGQJ33gk33givveY6GpHsF41GiUajaTlWQVqOUrWuwBTgSO/+QiACrAM6ANOBXkleF4vFKs0VOWHy\nZHj1VfjLX1xHklm7dkFxsW1/8gn07+82HpFcUlBQALU8t7tosHgRGO9tjwdecBBDKORTp7Jfo0aW\nBE86CR591HU0IhIXdEJ4Cngb6AmsBL4P3AmMxIadnuTdz0v5Muw0mUsugT/9Cf7+d9iUOOxARJwI\nug/hokr2nxLw+2aFfK0Q4kpLYdw4qxSmTYM2bVxHJJLf8vT7aTjkY6dyorvvhlGj4PjjbQTSihWu\nIxLJX0oIDuXjsNNEBQVw++02YW39etsWETfy/HTkVr43GcUVFMCFF8JDD8HTT8PSpfCTn8CWLa4j\nE8kvruch5LV87lROpnNnGDrU+haKi2HRIs1VEMkknY4cUoVwsHvvtUlrGzbAZ5/Bu++6jkgkfygh\nOKQK4WA9elg/QtOmcN11cP/9riMSyR86HTmkCqFqF18Mr7wCb77pOhKR/KCE4JCGnVatTRt48kn4\n3vfghRdsmY/t211HJZK71KnskIadVu/MM+GWW+Ccc6BXL2jZEl5/HRo3dh2ZSO7R6cghNRml5sor\nIRqFuXOhbVt44AHXEYnkJlUIDqlTOTWFhXDiibZ9xx129bUhQ+Doo6F+fVULIumi05FDqhBqrndv\nWyF1/Hjo2BEGDYJ9+1xHJZIblBAcUqdy7Zx9Nnz+OaxeDS1awEsvuY5IJDcoITikTuXaq1cPSkrg\nxz+GBx+EGTOsCUnVgkjtqQ/BIVUIdTd2LFx1FZxxhvUlRKNWPbRvD2PGuI5OJLsoITikCqHuiovh\n2WchFoOZM+EPf7A5CwC7d0ODBm7jE8kmOh05pE7l9DjtNDj9dFsY74UX4LHHYNgweOYZ15GJZBdV\nCA5p2Gl6nXaaNReVlsIRR9gM55ISSxYiUj2Xp6PlwKfALOB9h3E4owohvQoLLRkADB8Ov/kNXHaZ\nrYV06aXwi1+4jU8k7FxWCDEgAmx0GINT6lQO1ve/b0tdnH22fdZ79sDOndC1K/zsZ/acggKnIYqE\niusmo7z+c1SncvDGjIFPP7XO53nz4J13rON55kzrcP77311HKBIeriuE14H9wB+BRxzG4oSajDLj\n0EPtZyRit9at4eqroUkTmDLF9jVr5jBAkZBw+f30eGAAMBq4GhjuMBYn1KnsxrhxMG2a9StccAH8\n+teuIxIJB5cVwlrv55fA88AgYIb/CWVlZd9uRyIRIpFIhkLLDFUIbhQXW1UwfDjcdJNdpW3uXFsG\n4/XX7aI8Q4a4jlIkNdFolGg0mpZjuWrDbwwUAduAJsBrwC3ez7hYLBZzEFrmXHIJjB5t31TFnTlz\nbHLb22/DRRfBbbfBBx9YM9LKlXD44a4jFEldgY2UqNW53VWF0A6rCuIx/JUDk0FeUIUQDkceabe4\nhQuhW7eK+xddBH/8o0YkSe4L83/xnK8QLrgAzjvPfkq47N4Ny5dDly7WtFRUBJ07w//8jy27rb4f\nCatsrBAEDTsNs4YNoWdP2542DWbNsialHj2sUhg3DkaOtDkO9fRXJDlCpyOH1GSUHVq1gpNPhv/6\nL1izxuYzHHEE/Pa3cOqp8PXXNhv6V7/S8tuS3ZQQHNKw0+xTUmIzna+/3ia3fec7MGAAnHOOLaY3\ndqxVFHv3uo5UpOZU7DqkCiG7FRXBxIl2LYZmzezynhMn2ryGJUtg8GBo2hR++UvYvNlWYBUJM30/\ndUgVQm4YPtyu1tawoSWDd96Bp5+2DukVK2xZ7vPPt9VX27aFG26w1331Ffzud7a+kkgYqEJwSBVC\n7oovkxGL2e299+Caa+DJJ+EHP4ATTrDhrddfbwnhppv0BUHcU0JwSKud5r6CArsdd5xNdgP45z/h\nrLOsA/qJJywpfPghPP+8DUPu1Anuuefg/xuxmOZCSLD0fcQhDTvNT0OGwOLFdnW3yy+3zumRI23G\ndHk5PPKIrdL63e/CFVfAW2/B7bfb6qxpWqFAJClVCA6pySh/lZRYsxFAr152A3juOfjsM3j8cet7\nWLYMLrwQtm6Fhx6y+Q9dukD//jb09fDDYds2G/Z6yy02aU6ktsJcgOb8TOVhw+COO6xTUqQyO3fa\nzOmSEti1yyqGOXPscqHPP28X/rnwQuunePNNW9Zb8lddZiorITh03HHWVjx0qOtIJFtt22Z9C82a\nwZVXWqIYNMhmT/fvDyeeCB06WNLYtw8eftjmTQwbZqu7Su5RQshSgwfb2jiDB7uORHJBLAb33Qfr\n1llVsXYtvPGGTZIrL7c5ESNHWpPU2rX23LlzbdXdww6zJBL/k1PndfbSWkZZSp3Kkk4FBXDddQfu\n27QJ1q8OSryGAAAIB0lEQVS32+7dttQG2DWlzzvPrhx31FH2f3HECBsKW1AAN99scyr69YPJk625\nyp8kyssrRlBJ7gjzP2fOVwgDBsBjj8HAga4jkXwTi1kzUsOG9nPjRnj/fWje3CbTTZgAd94JL71k\n14o48kjYsMGSStOmVnUUF8PPf25JY88eGzXVr5/9f+7Tx57bpo3r3zT/qMkoSx11FEyaZLNcRcJq\n3z4b9dS7t41sWry4Yk7En/5kcyhWrrQlOpYutZna8+db5/bFF8Oxx9rzFy2yxDFlio2YGjHCksqa\nNZaYWre299u6FbZssdFUUnNKCFmqXz946qkDL84ikm3277cTfXGx3S8vtxVgt261k//LL9vJvqjI\nlhD/7W/hRz+y5qwePaw/o6jItgcNsol769bBo4/asaZOtcEXnTrBn/9so6uGDYNGjeynlh8/kBJC\nlurTx1bI7NvXdSQimRGvLPbutcrjww9txdhVq+CLL2DGDKsoCgrgwQctUTRubE1XLVrYaKqxY60K\n2bIFVq+uWDRw717r6zjtNEsW5eXWJ9KzpzWFvfIKjBpliWvwYKhf3+1nERQlhCzVsyf87/9WTEoS\nkeTiI6fat6+oRMCuajdjhiWOJk3g3Xdt5vfixXDIIZYcFiywTvWTT7ZFB/fvt2TQvLltb9tmVcaG\nDfY32bu3XUJ12zZLQJs3W79IvXo2VPzrr22k1tixtu+bbyyO2bOtGbhfP1efklFCyFJHHGHl9BFH\nuI5EJLdUte5TLGZNV/EKoWFDq1Y6drTmqAUL7ATftCls326VyZtvWsXx9tuWkLp1s+cWFlpC6NTJ\nOtBXroTRo62Zq08f6zcpLbW5ILGYVSznnmuLGQYlWxPCKOA+oAh4FJiY8HjOJ4Tu3e1iKqWlriMR\nkVTs32/Jo0EDO/k3bGhVROPG9vhdd1kn+bhxNlqrb1+7NsaqVVZN9OoFLVsG2ypQl4TgShHwOdAV\nqA/MBnonPCeWDaZPn17r1x56aCy2bFnaQqlUXWLMJMWZXoozvbIlTqDW36RdTYsahCWE5cBe4O/A\n2Y5iqZNoHZafzNT693WJMZMUZ3opzvTKljjrwlVC6ASs9N1f5e3LK1rtVETCxNUI3pRKmjPPDDqM\nulu0CD76qHav/fprjaEWkfBw1fEwBCjDOpYBbgTKObBj+XNA3a0iIjWzBDjcdRA1UQ8LuivQgOSd\nyiIikidGA4uwSuBGx7GIiIiIiEhYjQIWAp8BExzH8jiwHpjj29cKmAYsBl4DSnyP3YjFvRA4NUMx\nAnQBpgPzgLnANd7+MMXaCHgPax6cD9wRwhj9ioBZwBTvfhjjXA58isX5vrcvjHGWAM8CC7B/+8Eh\njLMn9jnGb1uwv6OwxRl/33nYeelvQMOQxllnqUxYy6ThwAAOTAh3Ab/wticAd3rbfbB462Pxf07m\nhvW2B+KLaDfFmuJ6hzBWbz4n9YB3gWEhjDHuZ8BfgRe9+2GMcxl2IvALY5yTgB942/WAFiGNM64Q\nWIt90QpbnF2BpVgSAPgHMD6EcabFccArvvs3eDeXunJgQlgItPO223v3wbKwv6J5BRtN5cILwCmE\nN9bGwAdAX8IZY2fgdWAEFRVCGONcBrRO2Be2OFtgJ7BEYYvT71RghrcdtjhbYV/4WmLJdQowMl1x\nhi1TZMOEtXZYMxLez/g/Qkcs3jhXsXfFqpr3CF+shdi3lfVUNHGFLUaA3wH/iQ2FjgtjnDEscX0I\nXOXtC1uc3YAvgSeAj4FHgCYhjNPve8BT3nbY4twI3AOsANYAm7GmorTEGbaEkG2r2VW3bkimf5+m\nwHPAtcC2JLG4jrUca9rqDJyAfQNPjMF1jGcAG7B25Mrm6YQhToDjseQ/Grgaa+JMjMN1nPWAgcCD\n3s9vOLjqD0OccQ2AM4FnKonDdZylwE+xL34dsb/5S5PEUas4w5YQVmPtdnFdODC7hcF6rCQD6ICd\nPODg2Dt7+zKlPpYMJmNNRhDeWLcA/wKOIXwxDgXOwppjngJOwj7TsMUJ1s4N9g38eWyNsLDFucq7\nfeDdfxZLDOsIV5xxo4GPsM8Uwvd5Hgu8DXwN7AP+iTW1h/XzrJMwTljrysGdyvE2uRs4uPOmAVYm\nLyFzM8ELgCexpg6/MMXahoqRD8XAv4GTQxZjohOp6EMIW5yNgWbedhNgJtb2HbY4wf6te3jbZV6M\nYYwTbKHN8b77YYvzKGwkYbH3fpOw6jBscaZNmCasPYW10+3B+ja+j3XqvE7y4V03YXEvBL6bwTiH\nYc0xs6kYNjcqZLEeibUhz8aGSv6ntz9MMSY6kYpRRmGLsxv2Wc7GThDxv5WwxQl2EvsA+AT7Rtsi\npHE2Ab6iItFCOOP8BRXDTidhrQNhjFNEREREREREREREREREREREREREREQkV7WmYv7GWmxm7Sxs\n+Y/fO4xLREQcuhlb+lok74RtLSORMIhP7Y9QsXRFGTYr9N/YhWnGAndjs66nYsuugK3PFMVWIH2F\nivVlREJPCUEkdd2wFVrPAv6CLTvcH9gJnI4tIfAAcC62CNkTwH87iVSkFupV/xQRwZYMngrsx9YO\nKgRe9R6bgy2C2AO76M/r3v4ibC0skayghCCSuj3ez3Jgr29/Ofa3VIAtOjY0w3GJpIWajERSk8qS\nwYuAQ6i4RGF9bPlhkayghCBysJjvZ7JtOPiqUzGsajgPmEjFUuTHBRemiIiIiIiIiIiIiIiIiIiI\niIiIiIiIiIiIiIiIiEP/H7DvjyvIf8HIAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "transfer.plot(response='time')" ] }, { "cell_type": "code", "execution_count": 49, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEZCAYAAACXRVJOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xm8XPP9x/HXzb4nyCIEUSQNQmgr1DaIFI1UedBa+kOj\nvx9F+Pm1lajlolVKS2tptRESJZZIU2ubCBMVmhBiySIRCRHJTWLJRq4s8/vj8x1z7twlk3vnnPOd\nM+/n4zGPOcss3+/Mvec93+/3LCAiIiIiIiIiIiIiIiIiIiIiIiIiIiIiTbAF+FoIr7sYOCaE1xVp\nkmZxF0BitRj4HFgbuP0xzgKViYy71eU+oBr7Lj4BpgD7RFMsKXcKhPKWAYYAHQO34SG8T/MQXjOr\nwt2SIgPchH0XOwEfAPfGWiIpGwoEqc85wIvAzdgv1feA4wLrOwP3AB8BHwLXk/t7OgeYBvweWAVc\nA2wPPAGsBmYAvwL+7R5/J3BL3vs/DlxaT9nS7vnTgPXA7sDXgcnAx8A84NTA408AZgNrXFn/L7Du\nJ8AC97x/AD3reL+BwDJqBs/3gTfcdDNgBPCuq+/DwHaBx/4IeN+tu6KeOtVlA/AoNVsIOwGPASuw\n7+TiwLqDgFexz3g58Du3vDfW/fUTYCn2nQU/g9bAbW7dUuBWoJVbl8I+s8uAKvfccwLPbeizHQLM\nAj7Fvqv+BdZbRGKwiPr7ss8BvgSGYRvC87GNRdbfgT8BbYFuwHTgvwPP3QhciG0s2wAPAQ+66X7Y\nL98X3OO/5V47u8Htim3ou9VTtjTW3dXPvX5nYAlwtpsfAKzEQgJsY36om+4MHOCmj3aPG4BtAP8I\nTA28T3AM4V1gUGDdo8Av3PQlwEvYxrol8GdXV4C9se6fw9x7/M59NkfXU7d7sXAFaA/cDzzn5psB\nM4ErgRZYEC4EBrv1LwNnuul2WJBBLhAewL6vfbFAyX7317nyd3W3aW4ZWCBsBCqxlt7x2HfT2a2v\n77M9AAuQb2Hf639hf2/ZoBERzyzGNlafBm7D3LpzsF/OWe2wjUp3oAf267VNYP3p5DZc52C/iLOa\nY+GyV2DZ9eRaCABzyG1wLwKebKDcz2MbqKwfkAuXrLuBq930+1hYdcp7zD3AjYH59q6cu7r5YCBc\n7x4P1p2zDtglUPbgBr6ne53mrgwPBta1w8YI6guE+4AvsO9iM7bB7+rWDaTm5wowEhjtpqdin0vX\nvMf0dnXpE1h2EzDKTS+kZutvMLbxBguEz6nZm1CFtUag/s/2T+RCJWsecATiLXUZlbcM8D2seyN7\nuyewfnlg+nN33wHYDfslvIxckPyZmr/olwSmu2G/aIPLPswry1jgLDd9FvbLGPe62QHvEfW8/m7Y\nxjIYbGdgwQVwCta1sRhrXRzslvek5gZ2PdZ1tDO1jQNOxn7hnoz9Us+WoTfWYsq+9xxgk3v/nnl1\n/dy9R30yWDfddu51q7Ff19l67pRXz5FYSIOFeR9gLtYt99281w5+Zh+Q6x7L/xw+cO+T9TEWKME6\ndHDT9X22u2HdR8Gy9qLuLjnxRIu4CyAlaQm2odqBmhuKoOBeNCuxDeQu5Fodu+Q9/m/AW8D+WFfP\nRLf8fHdr6PU/wH4dD67jcWD96idhv9gvBh7BWgEfYRvdrPZYnZZS2xxso3k8FjbBX/0fAOdiXTb5\nlmFdW1nt3Hs0JNt1tgQb5H8U+zW/BPvl3qee573ryga2oR6Pjd1k7Qq8E5j+yE1nP4e5dazbmvo+\n2w+AXwM3FPg64gG1EKQxe+gsAyZhg8Ydsb+jPai/O2AzMAHrzmiLbfB/RM2N+ofYxmUstiGr3oZy\nP4ltJM/CWi4tsb7rr7vpM7H+7c1YS2Oze944bEO+PzawegPwH2xjVpcHsYHuw7GNdNaf3XOzXU3d\ngKFuejw2uHoo1rq4job/7/K/j2exDf0F2DjNWmzsoi22Ed4X+KZ77FnkWmmrsc83GNhXuuftg3Xr\nPeyWj3PrsmMIV5NroTWkoc/2r1iQH+Tq1B5rsXSo/TIi4oNF1D4O4TG37mxq98tvJten3gm4C/vV\n+hnwGnBaA8/tim24V2MbthuxjV3QWdgG7MitlPt54Md5y/q411+B7c3zLLAfttF6BttTKvve3w48\n73+wDe7H2J5Nwa6SYH3BWjWbsb2lgiqA/8X6yNe41/tVYP1/UXMvo/doeFA5v+/9NKzV0hLrcnkQ\nC+VPsMHg7Gvdj/Xvr8VaW9lQ6o19rue511kG/Czw+q2BP2Ctgo+wPY6CexnlB+Qi955b+2y/g3Vd\nfepe92EUCGWvC/YraS7W7B6INWMnA/OxX5pdYiudxOUmau9ffzi1B02l6XpjgaAeAYndGHK/5lpg\nzcvfkttl73Jq7ukhydQX+8VegXUjrCT3Cxbs1+ZjWNeFFFdvFAjigc5Y8zjfPHJ7gOzo5iXZvokN\nKK/H/iYuD6zrh+3G+SLqUghDb6yrS4EgsRqA9Svei/Ux/xUbXPo08JiKvHkREYlB2L8YWgAHYoOP\nB2K/DkfkPaahE32JiEhEwj4O4UN3e8XNj8cOpFmOdRUtx/aaWJH/xD322COzcOHCkIsnIpI4C4E9\nG/PEsFsIy7HdErMH0gzCToT1BLZrIu5+Yv4TFy5cSCaTSeztmmuuib0Mqpvqp/ol74YdE9QoURyp\nfDF2Uq1WWHKdix1Q8wh2qP1icvuvi4hITKIIhDewo0bzDapjmYiIxES7ocUklUrFXYTQJLluoPqV\nuqTXryl8vtJUxvWHiYhIgSoqKqCR23a1EEREBFAgiIiIo0AQERFAgSAiIo4CQUREAAWCiIg4CgQR\nEQEUCCIi4igQRCRyFRWwalXcpZB8CgQRidSmTXb/+uvxlkNqUyCISKTWrrX7mTPjLYfUpkAQkUit\nXm33uv6VfxQIIhKpNWvsfvPmeMshtSkQRCRSCgR/KRBEJFLZQNiyJd5ySG0KBBGJVHYMQS0E/ygQ\nRCRSa9ZAx44KBB8pEEQkUmvWwHbbqcvIRwoEEYnUmjXQpYtaCD5SIIhIpDZsgA4dFAg+UiCISKQ2\nbID27dVl5CMFgohEasMGaNdOLQQfKRBEJFLZFoICwT8KBBGJVLaFoC4j/ygQRCRSaiH4S4EgIpHS\nGIK/WkTwHouBNcBmYCNwELA98DCwm1t/GvBZBGURkZipy8hfUbQQMkAKOAALA4ARwGSgDzDFzYtI\nGVCXkb+i6jKqyJsfCoxx02OAkyIqh4jETMch+CuqFsKzwKvAT9yyHkCVm65y8yJSBqqrNYbgqyjG\nEA4FlgHdsG6ieXnrM+4mImVAXUb+iiIQlrn7lcDfsXGEKmBHYDnQE1hR1xMrKyu/mk6lUqRSqRCL\nKSJRUJdRcaXTadLpdFFeK79vv9jaAc2BtUB7YBJwLTAI+Bi4CRtQ7kLtgeVMJqOGg0jSdOoETz8N\nF1wAb70Vd2mSp6KiAhq5bQ+7hdADaxVk3+sBLBReBR4BhpHb7VREyoC6jPwVdiAsAgbUsfwTrJUg\nImVk82bYtAnatFGXkY90pLKIRKa62sKgeXO1EHykQBCRyGzYoEDwmQJBRCKTDYRmzdRl5CMFgohE\nZsMGaN1aLQRfKRBEJDLBLiO1EPyjQBCRyAS7jNRC8I8CQUQio0FlvykQRCQy6jLymwJBRCKTPQ5B\nXUZ+UiCISGTUZeQ3BYKIREZdRn5TIIhIZLSXkd8UCCISGXUZ+U2BICKRUZeR3xQIIhKZ/HMZ6RpY\nflEgiEhksoFQUWE3BYJfFAgiEpnsye1AA8s+UiCISGSyLQTQwLKPFAgiEpn8QNDAsl8UCCISmWAg\nqMvIPwoEEYnM+vXQvr1Nq8vIPwoEEYnM+vXQoYNNq8vIPwoEEYnMunW5FoK6jPyjQBCRyKxbV7OF\noEDwiwJBRCKjLiO/KRBEJDLqMvKbAkFEIpPfQlAg+EWBICKRyGRqtxDUZeQXBYKIRKK6Glq2hBYt\nbF4tBP9EEQjNgdeBJ9z89sBkYD4wCegSQRlEJGbBPYxAg8o+iiIQLgHmANkT3Y7AAqEPMMXNi0jC\nBbuLQIPKPgo7EHoBJwCjgAq3bCgwxk2PAU4KuQwi4oHggDLoegg+CjsQbgV+DgQbhj2AKjdd5eZF\nJOGWLoUegf92BYJ/WoT42kOAFdj4Qaqex2TIdSXVUllZ+dV0KpUilarvZUTEd3PnQr9+uXkFQnGk\n02nS6XRRXqti6w9ptBuAHwGbgDZAJ2AC8C0sIJYDPYHnga/X8fxMRn8tIolx/vmw775w0UU2378/\nPPAA7LdfvOVKmoqKCmjktj3MLqMrgF2A3YEfAs9hAfE4cLZ7zNnAxBDLICKemD8f+vbNzauF4J8o\nj0PIfvU3Asdiu50e7eZFJOFWrYLu3XPzCgT/hDmGEDTV3QA+AQZF9L4i4olPPoHtt8/NKxD8oyOV\nRSQSCgT/KRBEJHRffGEHobVrl1umQPCPAkFEQvfpp7DDDhYCWQoE/ygQRCR0+d1FoEDwkQJBREKn\nQCgNCgQRCV1VFXTtWnOZAsE/CgQRCd3cufD1vPMRKBD8o0AQkdDNng17711zmQLBPwoEEQnd7Nmw\nzz41lykQ/KNAEJFQrV0LixcrEEqBAkFEQjVjBgwYAK1a1VyuQPCPAkFEQjVzJnzrW7WXKxD8o0AQ\nkVDNn197DyNQIPhIgSAioZo/H/r0qb1cgeAfBYKIhEqBUDoUCCISmjVrYN062Gmn2usUCP5RIIhI\naBYsgD33rHmW0ywFgn8UCCISmvq6i0CB4CMFgoiE5p13YK+96l6nQPCPAkFEQvPaa3DAAXWvUyD4\nR4EgIqHIZODll+Hgg+ter0DwjwJBREIxZ45dQ7lXr7rXKxD8o0AQkVA8/DCcckr96xUI/mkRdwFE\nJHmuugp+9SuYPr3+xygQ/FNIC2HvOpalilwOEUmI2bMtDKDuk9plKRD8U0ggPAJcDlQA7YDbgRvD\nLJSIlK6nn4aLL7ajlOs6IC1LgeCfQgJhILAL8DIwA1gGfDvMQolI6ZoxAw46CDp2bPhxCgT/FBII\nm4AvgLZAG+A9YEuYhRKR0jVjRsNdRVkKBP8UEggzgA3AN4HDgTOARwt4XhtgOjALmAP8xi3fHpgM\nzAcmAV22rcgi4qtFi6C6uv7TVQQpEPxTSCCcB1wFbMS6i4YCTxTwvA3AUcAAYD83fRgwAguEPsAU\nNy8iCTBlChx1VMNjB1kKBP8UEghVwK55t6kFvv7n7r4V0Bz4FAuUMW75GOCkQgsrIv7KZODuu+GH\nPyzs8QoE/xRyHMLTQPZrawPsDrwD7FPAc5sBrwF7AH8CZgM9sJDB3ffYhvKKiKdefBFWr4YTTyzs\n8QoE/xQSCPvmzR8IXFjg62/Buow6A//Cuo2CMuTCppbKysqvplOpFKlUqsC3FZEoTZkC118Pw4dD\nswLPf6BAKI50Ok06nS7KaxXQ01ent6kdFFtzFba30nnYgW3LgZ7A80Adl+Amk9Ffi0hJyI4ZrFgB\n3boV9pyTT4Yzz2z49Bay7Srsy2jUtr2QFsL/BaabYS2EpQU8ryu2y+pn2C6rxwLXAo8DZwM3ufuJ\n21BeEfFMJgMtW8JddxUeBqAWgo8KCYSO5Lp1NgFPAo8V8Lye2KBxM3e7H9ur6HXs6OdhwGLgtG0q\nsYh45fzzYeNGOO+8bXueAsE/hQRCZSNf+y2sNZHvE2BQI19TRIpk0SLrshk0CK67btufv2EDfPe7\n8NxzsN122/58BYJ/GgqEho41yGC7j4pIifra1+y+urpxgXD77fDqq/CXv8CwYdv+fAWCfxoKhN81\nsE5fo0hCtG+/bY+fOhUGDLCT2D30EBx/fOPet6ICtugkOF5pKBAWAe9HVRARiUfbtoU/9osvIJWC\niy6y1sFhhzX+fZs1UwvBNw3tMRzc+6eQQWQRKRHBX+Zr1uSWffhhw8+bNMnu77gDRo7c+hlNG6Iu\nI/8UegnNr4VaChGJ1GefQadO8I9/QFUV/OIXMHkynJa3z18mY8cWZE2fDiecYNOXXda0MigQ/KNr\nKouUoZUroUcPOxHdokVw882QTsOqVbBpEwwcCF9+CS+9ZI+rqrIWxLRpcMEFtodRmzZNK4MCwT8N\njSHsB6x1020D02CDyp3CKpSIhGv5cjuILNjlM306fPopLFli1zRYvBg++MDWjR9vQfHZZ3DEEdC6\nddPLoEDwT0OB0DyyUohIpJ5/Hg45pOay6dOtVbBwoc1fc43tRdSsGcyda62D226zrqZiUCD4R11G\nImXoySdhqDuSaOxY+OY34fPPrRXw5pu2/KGH7H7wYDt53QcfWOugWBQI/lEgiJSZ6mqYMyd3mcsf\n/QjuvDO3fuZM2Gmn3Py++8K8edCvHzQvYr+BAsE/hZy6QkQS5O23Yc89ax5/0L9/bnrePDjuOBg9\n2pZffDF0726BUEwKBP+ohSBSZh57DA4+uOayYDgsWGDnNwJrOey6K/z85zBkSHHLoUDwjwJBpIyM\nHQt/+xtce23tdevW2UZ/7VrYfXc46CD4el1XKikSBYJ/1GUkUiYyGbjlFhgzBnr2rL2+fXvYZReb\n7tbN9joKkwLBP2ohiJSJqVPtXERHHln/Y/r2tfttudBNYykQ/KMWgkgZWL4czjgD/vCHhq95vNtu\ndt+UcxQVSoHgH7UQRMrAtGl2rMGppzb8uOyeRBWNvdr6NlAg+EeBIFIG3nwT9t9/64/r2ze6jbQC\nwT8KBJEy8OqrhQVClBQI/lEgiCTcggV2srrBg+MuSU0KBP8oEEQSbtQo+PGPi3dSumJRIPhHexmJ\nJNz48TBhQtylqE2B4B+1EEQSbOVK+Pjjmucq8oUCwT8KBJEEe+kl+MY3Gj72IC4KBP94+GciIsXw\nhz/YNZLPOy/uktRNgeAfjSGIJNS998KkSQ2fqiJOCgT/qIUgkkBLl9oVzg49NO6S1E+B4B8FgkjC\nPPMM9OoFRx8NLTzuA1Ag+CfsQNgFeB6YDbwNDHfLtwcmA/OBSUCXkMshUjYuvRRatcpdM9lXCgT/\nhP37YSPwv8AsoAMwEwuCc939b4HLgRHuJiJN8NFHsGqVnebaxz2LghQI/gn7T2Y5FgYA64C5wM7A\nUGCMWz4GOCnkcoiUhalT4Ygj/A8DUCD4KMo/m97AAcB0oAdQ5ZZXuXkRaaJ0GlKpuEtRGAWCf6Ia\ncuoAPAZcAqzNW5dxt1oqKyu/mk6lUqRK5S9dJAabN9tupj/9adwlKYwCoTjS6TTpdLoorxXBZTBo\nCTwJPAPc5pbNA1JYl1JPbOA5/3LemYz+WkQKduWV8MIL1m0UxQVummrECOjcGUaOjLskyVJhX36j\n/gLC7jKqAO4B5pALA4DHgbPd9NnAxJDLIZJo77xjZzUdN640wgDUQvBR2F1GhwJnAW8Cr7tlI4Eb\ngUeAYcBi4LSQyyGSaBMmwCmnwM47x12SwikQ/BN2ILxI/a2QQSG/t0jZmDLFjj8oJQoE/5TAzmki\n0pA33rBA8Pk0FXVRIPhHgSBSwpYvh6OOguHDYbvt4i7NtlEg+EeBIFLC/vEPOOEEO9V1qVEg+EeB\nIFKiVq+2IDjmmLhL0jgKBP8oEERK1K9/bRvVIUPiLknjKBD84/HJcUWkPtXVcP/98Nxz0K1b3KVp\nnIoK2LIl7lJIkFoIIiUmk7HWwf77Q79+cZem8dRC8I9aCCIlZtQoeOwxeOqpuEvSNAoE/ygQRErI\nli1w443w4IPQu3fcpWkaBYJ/1GUkUkJeeAHatYODDoq7JE2nQPCPAkGkRKxaZWcIPffc0jmBXUMU\nCP5RIIiUiLFjoXt3uPDCuEtSHAoE/2gMQaQEvPAC/OxndlbT1q3jLk1xKBD8oxaCiOcyGTtX0U9/\nCscfH3dpikeB4B+1EEQ898orsH493H57MsYOshQI/lELQcRzN90Ew4YlKwxAgeAjtRBEPDZxIrz9\nNjzwQNwlKT4Fgn/UQhDx2LXXwq23Qps2cZek+BQI/lEgiHjqlFPsAjjHHRd3ScKhQPCPuoxEPLR0\nKTz7LCxZAs0S+rNNgeCfhP6piZS2e+6B00+HTp3iLkl4FAj+UQtBxDMLFsAdd0A6HXdJwqVA8I9a\nCCKeufpquOwy2HvvuEsSLgWCf9RCEPHIrbfCiy/C3XfHXZLwKRD8o0AQ8cSCBfCb38BrryV77CBL\ngeAfdRmJeGD1avjBD2DIEOjVK+7SREOB4B8FgogHHnwQXn8dzjgj7pJER4HgHwWCSMyqq2H0aBg/\nHgYNirs00VEg+CfsQBgNVAFvBZZtD0wG5gOTgC4hl0HEa3fdBV27wve+F3dJoqVA8E/YgXAvkH/g\n/QgsEPoAU9y8SFkaNw6uvNIujdmizHbxUCD4J+xA+Dfwad6yocAYNz0GOCnkMoh4af58u/DNxIlw\n5JFxlyZ6CgT/xDGG0APrRsLd94ihDCKxymTgiivgoovg2GPjLk08FAj+ibuRmnE3kbKRydhFb957\nD+6/P+7SxEeB4J84AqEK2BFYDvQEVtT3wMrKyq+mU6kUqVQq5KKJhOv996F3b5tesADato21OLFS\nIBRHOp0mXaQTX0VxUb7ewBNAfzf/W+Bj4CZsQLkLdQ8sZzL6a5GEufZa+P3vYeNG+PzzuEsTr9Gj\n4d//hnvvjbskyVJh11pt1LY97DGEccBLQF9gCXAucCNwLLbb6dFuXiTx/vMfeOopePhhWL8+7tLE\nTy0E/4TdZXR6PcvL6PAbERg1Cq66CvbfHw47zDaG5U6B4J+4B5VFEm/WLDvWYPJk6N9/648vFwoE\n/+jUFSIh2rQJ/vhHu/qZwqAmBYJ/1EIQCdH999ug6dSpcZfEPwoE/6iFIBKSTAbGjrVB5COOiLs0\n/lEg+EeBIBKCN96AH/7Q9iYaOjTu0vhJgeAfBYJIkX32mZ2OIpOBJ5+ENm3iLpGfFAj+USCIFNGj\nj8Iee8Ahh8Ajj0D37nGXyF8KBP9oUFmkSN58Ey680K56dnp9R+DIVxQI/lEgiDTR+vXw0kswciTc\ncAOcd17cJSoNCgT/qMtIpIkuuggGD4bNm2HYsLhLUzoUCP5RC0Gkkd591y59OXEifPQR7LCDTkmx\nLRQI/lEgiDTCli2w117Qvr21Dnr2jLtEpUeB4B91GYlsow8/tBBo3twGkK+5Ju4SlSYFgn8UCCLb\n4JVXYMAAu/3rX/CXv9gZTGXbKRD8oy4jka1YswZWrrTjCy65BG67Dc46K+5SlT4Fgn/UQhDZiiuu\ngD33hDvugBUr7JQU0nQKBP+ohSDSgDVrYMIEOPVUuPhiu+pZC/3XFIUCwT9qIYjUYfNmSKfhwAPh\n+9+HceNs19KBA+MuWXIoEPyj3zoiAdXVUFkJ991nG6vRo+H4423jpV1Li0uB4B+1EESwFsHw4XZw\n2dy5NnC8bBmccIIONguLAsE/aiFIWXv3XXjxRXj5ZQuCxYvt6GMJnwLBPwoEKVtr18KJJ9r0AQfA\n449Dly7xlqmcKBD8o0CQsrJ6NVx3nbUCRo+G73zHdieV6CkQ/KNAkLKxahVceqmdemKPPSwIvvOd\nuEtVvlq1skF88YcCQRJt2jS7hOUtt9gFbObMgeefh1Qq7pJJ5852nIf4Q4EgifPll/DEE3ah+5tv\ntvkrr4TTToMhQ6Bly7hLKGCBsHp13KWQIAWCJMbChfDPf8KoUdChA+y4ox1c1qkT9OsXd+kknwLB\nPz7vYZ3JaMRJtmLZMjuI7I47rCVw4olw9NFw5pk6fsB31dXQsaPd67sqngr7MBv1icbZQjgOuA1o\nDowCboqxLFIi1q2zM4/ed5+NBzzzjB1BPH48dO9ug8VSGlq3hmbNYMMGaNs27tIIxHekcnPgDiwU\n9gZOB8qqUZ9Op+MuQmiKWbf162HGDLj6ajvLaN++cPDB8NprcPjhtsfQggVwyCHRhUGSvzuItn5x\nDCwn/ftrirgC4SDgXWAxsBF4CPheTGWJRZL/KJtat/nzYdIk+OUvLQAGDoTrr4f99rMWQVWVDRoP\nHx7PgWRJ/u4g2vp16hT9OELSv7+miKvLaGdgSWD+Q0DnkSxzy5fbWUUrK6F/f2sB3HknzJoFw4ZB\nr15xl1CKTQPLfokrEAoaLc6eViCJ3nkHZs6MuxThaEzdMhnrBho82K4/cMwxuXXfK6u2Y3nZbju7\nCt0OO0T3no35+zz8cPjFL8Ipj0/iGts/GKjExhAARgJbqDmw/C6gIUIRkW2zENgz7kJsixZYoXsD\nrYBZlNmgsoiI5BwPvIO1BEbGXBYREREREfHVccA8YAFwecxlaazRQBXwVmDZ9sBkYD4wCQjuMDkS\nq+88YHBEZWyKXYDngdnA28BwtzwJdWwDTMe6MecAv3HLk1C3oObA68ATbj5J9VsMvInVb4ZblqT6\ndQHGA3Oxv9GBJKt+X2mOdSH1BlpSumMLhwMHUDMQfgtk91O4HLjRTe+N1bMlVu938f/SpjsCA9x0\nB6zrrx/JqWM7d98C+A9wGMmpW9ZlwAPA424+SfVbhG0gg5JUvzHAj910C6AzyarfVw4B/hmYH+Fu\npag3NQNhHtDDTe/o5sHSO9gS+ie2F1YpmQgMInl1bAe8AuxDsurWC3gWOIpcCyFJ9VsE5O/ImpT6\ndQbeq2N5UernW1LUdcDazjGVpdh6YN1IuPvsl7cTVs+sUqtzb6w1NJ3k1LEZ9quqilzXWFLqBnAr\n8HNsV++sJNUvgwXeq8BP3LKk1G93YCVwL/Aa8FegPUWqn2+BUC6nN83QcF1L5XPoADwGXAKszVtX\nynXcgnWJ9QKOwH5JB5Vy3YYAK7D+9fqOQyrl+gEciv1IOR64EOvCDSrl+rUADgTucvfrqd2L0uj6\n+RYIS7EBy6xdqJlupawKa8oB9MT+KaF2nXu5Zb5riYXB/ViXESSvjquBp4BvkJy6fRsYinWrjAOO\nxr7DpNQPYJm7Xwn8HTt3WlLq96G7veLmx2PBsJxk1K+GJB2w1pvag8rZvrwR1B70aYU1Bxfi93Uq\nwMo3Fut6CEpCHbuS20OjLfACcAzJqFu+I8mNISSlfu2Ajm66PTAN27MmKfUD+5vs46YrsbolqX41\nJOGAtXGHe0HkAAACCklEQVTAR8CX2JjIudheD89S925hV2D1nQeUwmXfD8O6VWZhXQ+vY7sLJ6GO\n/bG+2VnYros/d8uTULd8R5Lbyygp9dsd++5mYbtEZ7chSakfwP5YC+ENYAI20Jyk+omIiIiIiIiI\niIiIiIiIiIiIiIiIiN9K6gAFkSLYjB1fkDUOO6hHRETKTP45l4qhRQivKRI5385lJBKXxdhpAGZi\nLYi+bnl77IJH07EjmIe65edgR/lOwS5M0hZ4BDsz6gTsOgrfwI5SD57i4yfA78OqhIiIFG4TudNt\nvA6c6pYvws6MCXABdlphgBuAM910F+y0Ku2wQFhC7hQBPwP+5Kb3ATZiJx1rj502oLlbN82tF/GO\nmrpSbr7ATo1clwnu/jXgZDc9GDgR2+ADtAZ2xU4hPBn4zC0/FLjNTc8mN06xHnjOvcY87Cyxs5ta\nCZEwKBBEcqrd/WZq/m+cjF2TNmggtrEPqm8njVHAL7Fr4I5uYhlFQqMxBJGG/QsYHpjPti7yN/7T\ngNPc9N7YWVOzZmDnoT8D26tJxEsKBCk3bak5hnBDHY8JXnHqeqyb503sdMrX1vEYsCtYdcO6g653\n96sD6x8BXsxbJiIiCdQMG18A2AO7EHqw2+kJal+KU0REEqgjdtGSWdiFS7IXIsnumfRwTOUSERER\nERERERERERERERERERERERERKab/B5WvRS57HCkNAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "transfer.plot(response='energy')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "By enabling `save=True` parameter, the plots can be also saved." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# IO" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "TransferFunction can be saved in pickle format and retrieved later." ] }, { "cell_type": "code", "execution_count": 50, "metadata": { "collapsed": true }, "outputs": [], "source": [ "transfer.write('transfer.pickle')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Saved files can be read using static `read()` method." ] }, { "cell_type": "code", "execution_count": 51, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([ 0., 0., 0., 0., 0., 0., 0., 0., 0.])" ] }, "execution_count": 51, "metadata": {}, "output_type": "execute_result" } ], "source": [ "transfer_new = TransferFunction.read('transfer.pickle')\n", "transfer_new.time[1:10]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Artificial Responses" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For quick testing, two helper impulse response models are provided." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1- Simple IR" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "simple_ir() allows to define an impulse response of constant height. It takes in time resolution starting time, width and intensity as arguments." ] }, { "cell_type": "code", "execution_count": 52, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 52, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEACAYAAAC08h1NAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEW1JREFUeJzt3W+MVNd9h/FnvAtpUmRTpAjXsOpaQBSo2gYrosRNlGlD\nW4ocaCMlGMVxa0c2L0rttFFCcKV6LfVF0sqKi5AxTUlCaydUxZaFJWzaKJ4qUVxsYoxtYClLSwNY\n2JZTpzGKFXZ3+uLcZYfxLnN3ZnaHO7/nI43m/jn3zjns8L1nzz0zC5IkSZIkSZIkSZIkSZIkSSqA\n1cAgcALYPMH+9wPPAG8Dn6/Z3gc8DRwBXgbumt5qSpLaoQcYAvqBWcALwNK6Mu8FPgj8NZcG/7XA\nB7LlOcDxCY6VJM2gq3KUWUEK/lPABWA3sK6uzOvAwWx/rXOkCwXAW8Ax4Lom6ypJaoM8wb8AOF2z\nfibbNlX9wHLgQBPHSpLaJE/wV9vwOnOAPcDdpJ6/JKlDenOUOUu6STumj9Trz2sW8CjwMPB4/c5F\nixZVT548OYXTSZKAk8DiZg7M0+M/CCwhDdXMBtYDeycpW5pgfSdwFHhgogNOnjxJtVrt2se9997b\n8TrYPtsXsX3d3LZqtQqwKEd+TyhPj38Y2ATsJ83w2Um6Sbsx27+DNHvnOeBqYJQ0pLOMNKPnFuBF\n4FBWfgvwVLMVliS1Jk/wAzyZPWrtqFk+x6XDQWO+T77fKiRJM8RQnmblcrnTVZhWtq/Yurl93dy2\nVtWPyXdCNRuvkiTlVCqVoMkMt8cvScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEv\nScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEvScEY/JIUjMEvScEY\n/JIUjMEvScEY/JIUjMEvScEY/JIUTJ7gXw0MAieAzRPsfz/wDPA28PkpHitJmmGlBvt7gOPAKuAs\n8BywAThWU+a9wK8Afwj8L3D/FI4FqFar1eZbIEkBlUolaJzhE2rU418BDAGngAvAbmBdXZnXgYPZ\n/qkeK0maYY2CfwFwumb9TLYtj1aOlSRNk94G+1sZg3H8RrqCffaz8Nprna6FpqKvDx58sPXzNAr+\ns0Bf7euSeu555D52YGDg4nK5XKZcLud8CUnN+uY34bHH4Crn9hXCSy9VOH68Qk1cNq3RjYFe0g3a\njwGvAM8y8Q1agAHgp4zf3M17rDd3pRk2Ogo9Pem51NTtQXVaKzd3G/X4h4FNwH7SLJ2dpODemO3f\nAVxLmrFzNTAK3A0sA96a5FhJHTYyAr29hn5UV8KP3R6/NMN+9jOYNy89q5imczqnpC40PJx6/IrJ\n4JcCGhlJY/yKyeCXArLHH5vBLwVk8Mdm8EsBGfyxGfxSQAZ/bAa/FJDBH5vBLwU0POysnsgMfimg\nsU/uKiaDXwrIoZ7YDH4pIIM/NoNfCsjgj83glwIy+GMz+KWAnNUTm8EvBeSsntgMfikgh3piM/il\ngAz+2Ax+KSCDPzaDXwrI4I/N4JcCclZPbAa/FJCzemIz+KWAHOqJzeCXAjL4YzP4pYAM/tgMfikg\ngz82g18KyFk9sRn8UkDO6onN4JcCcqgnNoNfCsjgjy1P8K8GBoETwOZJymzN9h8Gltds3wIcAV4C\nvgW8q+maSmobgz+2RsHfA2wjhf8yYAOwtK7MGmAxsAS4E9iebe8H7gBuAH4tO9fN7ai0pNYY/LE1\nCv4VwBBwCrgA7AbW1ZVZC+zKlg8Ac4H5wP9lx7wH6M2ez7aj0pJaY/DH1ij4FwCna9bPZNvylPkx\ncD/wI+AV4E3gO61UVlJ7jIw4nTOyRtf8as7zlCbYtgj4HGnI5yfAvwCfBh6pLzgwMHBxuVwuUy6X\nc76spGbY4y+eSqVCpVJpy7ka/ejPAn01632kHv3lyizMtpWBHwBvZNsfA26kQfBLmn4Gf/HUd4rv\nu+++ps/VaKjnIOmmbT8wG1gP7K0rsxe4NVteSRrSeRU4nq2/m/QbwSrgaNM1ldQ2Bn9sjX70w8Am\nYD9pVs5O4BiwMdu/A9hHmtkzBJwHbsv2vQD8I+niMQo8D/x9G+suqUkGf2wTjc3PtGq1mvdWgqR2\nuOMOWLEiPauYSqUSNJnhfnJXCshZPbEZ/FJADvXEZvBLARn8sRn8UkAGf2wGvxSQwR+bwS8FZPDH\nZvBLAfmnF2Mz+KWA/NOLsRn8UkAO9cRm8EsBGfyxGfxSQAZ/bAa/FJDBH5vBLwXkrJ7YDH4pIGf1\nxGbwSwE51BObwS8FZPDHZvBLARn8sRn8UkAGf2wGvxSQs3piM/ilgJzVE5vBLwXkUE9sBr8UkMEf\nm8EvBWTwx2bwSwEZ/LEZ/FJAzuqJzeCXAnJWT2wGvxTM6Gh6vsr//WH5o5eCcXxfeYJ/NTAInAA2\nT1Jma7b/MLC8ZvtcYA9wDDgKrGy6ppLawuBXo+DvAbaRwn8ZsAFYWldmDbAYWALcCWyv2fd3wL7s\nmF8nXQAkdZDBr0bBvwIYAk4BF4DdwLq6MmuBXdnyAVIvfz5wDfAR4OvZvmHgJy3XWFJLDH41Cv4F\nwOma9TPZtkZlFgLXA68D3wCeB74GvKeVykpq3ciIUzmja3Tdr+Y8T2mC43qBG4BNwHPAA8CXgL+q\nP3hgYODicrlcplwu53xZSVNlj7+YKpUKlUqlLeeqD+x6K4EB0hg/wBZgFPhKTZmHgAppGAjSjeCP\nZud+htTzB/gwKfhvqnuNarWa9/oiqVWnT8ONN6ZnFVepVILGGT6hRkM9B0k3bfuB2cB6YG9dmb3A\nrdnySuBN4FXgHGkI6H3ZvlXAkWYqKal97PGr0Y9/mDRUs580w2cnaWbOxmz/DtKsnTWkm8Dngdtq\njv8z4BHSReNk3T5JHWDwq6lfE9rMoR5pBh07Bp/4RHpWcU3nUI+kLuMXtMngl4LxC9pk8EvBOMYv\ng18KxuCXwS8FY/DL4JeCMfhl8EvBOKtHBr8UjLN6ZPBLwTjUI4NfCsbgl8EvBWPwy+CXgjH4ZfBL\nwTirRwa/FIyzemTwS8E41CODXwrG4JfBLwVj8Mvgl4Ix+GXwS8E4q0cGvxSMs3pk8EvBONQjg18K\nxuCXwS8FY/DL4JeCMfhl8EvBGPwy+KVgRkaczhmdwS8FY49fBr8UjMGvPMG/GhgETgCbJymzNdt/\nGFhet68HOAQ80WQdJbWRwa9Gwd8DbCOF/zJgA7C0rswaYDGwBLgT2F63/27gKFBttbKSWmfwq1Hw\nrwCGgFPABWA3sK6uzFpgV7Z8AJgLzM/WF5IuDP8AlFqvrqRWGfxqFPwLgNM162eybXnLfBX4AjDa\nQh0ltZGzetToup93eKa+N18CbgJeI43vly938MDAwMXlcrlMuXzZ4pJaYI+/mCqVCpVKpS3najT8\nshIYII3xA2wh9d6/UlPmIaBCGgaCdCO4DNwFfAYYBn4BuBp4FLi17jWq1arD/9JM+eQn4VOfSs8q\nrlKpBE0OoTca6jlIumnbD8wG1gN768rsZTzMVwJvAueAe4A+4HrgZuC7vDP0Jc0we/xq9OMfBjYB\n+0kzfHYCx4CN2f4dwD7SDdwh4Dxw2yTnslsvXQEMfuX58T+ZPWrtqFvf1OAc/549JHWYwS8/uSsF\n459elMEvBeOfXpTBLwXjUI8MfikYg18GvxSMwS+DXwrG4JfBLwXjrB4Z/FIwzuqRwS8F41CPDH4p\nGINfBr8UjMEvg18KxuCXwS8F46weGfxSMM7qkcEvBeNQjwx+KRiDXwa/FIzBL4NfCsbgl8EvBVKt\nppu7zuqJzeCXAhkdhauuglKp0zVRJxn8UiAO8wgMfikUg19g8EuhGPwCg18KxeAXGPxSKAa/wOCX\nQnEqp8Dgl0Kxxy8w+KVQDH5B/uBfDQwCJ4DNk5TZmu0/DCzPtvUBTwNHgJeBu5quqaSWGfyCfMHf\nA2wjhf8yYAOwtK7MGmAxsAS4E9iebb8A/Dnwq8BK4E8nOFbSDDH4BfmCfwUwBJwiBfluYF1dmbXA\nrmz5ADAXmA+cA17Itr8FHAOua6nGkppm8AvyBf8C4HTN+plsW6MyC+vK9JOGgA5MrYqS2sVZPYJ8\nwV/Nea76r32qPW4OsAe4m9Tzl9QB9vgFkOctcJZ0k3ZMH6lHf7kyC7NtALOAR4GHgccneoGBgYGL\ny+VymXK5nKNakqbK4C+uSqVCpVJpy7nyfDlrL3Ac+BjwCvAs6QbvsZoya4BN2fNK4IHsuUQa+3+D\ndJN3ItVqNe8vFZJa8b3vwT33pGcVWyl9t3ZTX7Cd59o/TAr1/aQZPjtJob8x278D2EcK/SHgPHBb\ntu+3gFuAF4FD2bYtwFPNVFZSa+zxC/IFP8CT2aPWjrr1TRMc9338kJh0xTD4BYayFMrwsLN6ZPBL\noYyM2OOXwS+F4lCPwOCXQjH4BQa/FIrBLzD4pVAMfoHBL4XirB6BwS+F4qwegcEvheJQj8Dgl0Ix\n+AUGvxSKwS8w+KVQDH6BwS+F4qwegcEvheKsHoHBL4XiUI/A4JdCMfgFBr8UisEvMPilUAx+gcEv\nhWLwCwx+KZSREadzyuCXQrHHLzD4pVAMfoHBL4Vi8AsMfikUg19g8EuhGPwCg18KxVk9AoNfCsUe\nv8Dgl0Ix+AX5gn81MAicADZPUmZrtv8wsHyKx0qaIQa/oHHw9wDbSAG+DNgALK0rswZYDCwB7gS2\nT+HYrlepVDpdhWll+4qlPvi7rX21urltrWoU/CuAIeAUcAHYDayrK7MW2JUtHwDmAtfmPLbrdfub\nz/YVi8EvaBz8C4DTNetnsm15ylyX41hJM8g/vSiARqN91ZznKbVSiY9/vJWjr2zHj8MPf9jpWkwf\n21csL74Is2Z1uha60q0EnqpZ38I7b9I+BNxcsz4IzM95LKThoKoPHz58+JjSY4hp0gucBPqB2cAL\nTHxzd1+2vBL4jykcK0m6Av0BcJx0ddmSbduYPcZsy/YfBm5ocKwkSZKkKLrpA159wNPAEeBl4K5s\n+zzg34D/BP6VNN21yHqAQ8AT2Xo3tW8usAc4BhwFfpPuat8W0vvzJeBbwLsodvu+DrxKas+Yy7Vn\nCylrBoHfm6E6Nmuitv0t6b15GHgMuKZmX2Ha1kMaAuoHZlH8ewDXAh/IlueQhriWAn8DfDHbvhn4\n8sxXra3+AngE2Jutd1P7dgG3Z8u9pP9Y3dK+fuC/SGEP8M/AH1Ps9n2E9E0BteE4WXuWkTJmFunf\nYogr+ytrJmrb7zJe5y9T0LZ9iEtn/Xwpe3SLx4FVjM9ygnRxGOxYjVq3EPgO8NuM9/i7pX3XkIKx\nXre0bx6pM/JLpIvaE6QgKXr7+rk0HCdrT/2swqdIk1GuZP1c2rZafwQ8nC1PuW2dvCrk+XBYUfWT\nrtYHSG/CV7PtrzL+piyirwJfAEZrtnVL+64HXge+ATwPfA34RbqnfT8G7gd+BLwCvEkaEumW9o2Z\nrD3XkTJmTNHz5nbGZ1NOuW2dDP5qB197Os0BHgXuBn5at29s/m0R3QS8Rhrfn+wDe0VuXy9pRtqD\n2fN53vkbaJHbtwj4HKlTch3pfXpLXZkit28ijdpT1Lb+JfBz0n2ayVy2bZ0M/rOkG6Jj+rj0qlVE\ns0ih/0+koR5IvY5rs+VfJoVnEd1I+l6m/wa+DfwOqZ3d0r4z2eO5bH0P6QJwju5o3weBHwBvAMOk\nm4MfonvaN2ay92N93izMthXNn5A+O/Xpmm1Tblsng/8g6Rs9+0kf8FrP+A3DIioBO0mzQR6o2b6X\ndBON7Plxiuke0pvretIntb8LfIbuad850tDj+7L1VaQZME/QHe0bJI37vpv0Xl1Feq92S/vGTPZ+\n3Et6384mvYeXAM/OeO1as5o01LoOeLtme+Ha1k0f8Powaez7BdJwyCHSD2oe6YZoEafLTeajjF+k\nu6l9v0Hq8ddOl+um9n2R8emcu0i/oRa5fd8m3a/4OemifRuXb889pKwZBH5/Rms6dfVtu500XfN/\nGM+XB2vKF6ltkiRJkiRJkiRJkiRJkiRJkiRJktS9/h+lSqkVoCGKbwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "s_ir = simple_ir(dt=0.125, start=10, width=5, intensity=0.1)\n", "plt.plot(s_ir)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 2- Relativistic IR" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A more realistic impulse response mimicking black hole dynamics can be created using relativistic_ir(). Its arguments are: time_resolution, primary peak time, secondary peak time, end time, primary peak value, secondary peak value, rise slope and decay slope. These paramaters are set to appropriate values by default." ] }, { "cell_type": "code", "execution_count": 53, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 53, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHiBJREFUeJzt3XmUVOWZx/Fv0Q0Csq+yaBpZBCIEJMEmQSgWFTTAxNhB\nYtSgjrjOzMlMgphRO05iNCee5KhxCUZjjIIRSYRIQEHKBVFEkUAABRoUcGFRUNka6J4/nmq7qG6o\nqu5bde+t9/c5557ablc/tNav337ue98LIiIiIiIiIiIiIiIiIiIiIiIiIiKh8DDwMbAqxX7fAA4D\nF2S9IhERqVWDNPZ5BBiTYp8C4E5gPhCpb1EiIlI36YT6y8CnKfa5AZgF7Kh3RSIiUmfphHoqXYAJ\nwP3xx5UevKeIiNSBF6H+W+BGLMwjqP0iIuKbQg/eYxAwM36/HTAWOATMSdype/fulRs3bvTg24mI\nOGUj0CPdnb0YqZ8KdItvs4BrSAp0gI0bN1JZWRn47dZbb/W9Blfq7Nq1ktLSYNcYpE11ulkn0D2T\nQE5npD4DGI6NwrcAtwIN4689mMk3E6ly4ABs2wafpjoELyIZSSfUJ2XwfpPrWoi45b33oLJSoS7i\nNS/aL3klGo36XUJawl5nWRlEIrB7d27rqU3Yf5ZBozr9lcuZKpXx/pAIv/sd3HsvtG8PL73kdzUi\nwRWJRCCDrNZIXXxRVgaDBqn9IuI1hbr4YtMmOOMMhbqI1xTq4guN1EWyQz11ybnKSmjZ0mbAdOgA\ne/dCo0Z+VyUSTOqpS+Dt2gWFhdC6NbRqFYwZMCL5QqEuOVdWBt262f3WrdWCEfGSQl1yrqwMTj3V\n7rdqpVAX8ZJCXXJu06bqUNdIXcRbCnXJucSReuvW6qmLeEmhLjmXHOoaqYt4R6EuOadQF8kehbrk\n1KFD8MEHcMop9lgHSkW8pVCXnHr/fejcGRrGV+TXSF3EWwp1yalNm6rnqIMOlIp4TaEuOZXYTweN\n1EW8plCXnFKoi2SXQl1yKjnUdaBUxFsKdckpjdRFsktL70pOtWkD774L7drZ44oKmwlTXg4FBf7W\nJhJEWnpXAmv3bjh8GNq2rX6uQQNo0QL27PGvLpF8olCXnKmazhhJGnOory7inXRC/WHgY2DVMV6/\nGFgJ/BNYAvT3pjTJN8n99Crqq4t4J51QfwQYc5zXy4BhWJj/H/B7D+qSPKRQF8m+dEL9ZeB4H7ml\nQFVH9HWga32Lkvx0vFDXWaUi3vC6p34FMM/j95Q8oZG6SPYVevheI4DLgW8da4fS0tIv70ejUaLR\nqIffXoKsshJWrYLevWu+pgOlItVisRixWKzOX5/u3MciYC7Q7xiv9wdmY733DcfYR/PUHbZ2LZx3\nno3Wk2e/3H47fPYZ3HGHP7WJBJkf89RPwQL9Bxw70MVxCxfCqFE1Ax3UfhHxUjrtlxnAcKAdsAW4\nFYivhs2DwC1Aa+D++HOHgMHelilht2gRTJxY+2s6UCriHS0TIFl3+DC0bw/r1kHHjjVfX7AA7roL\nnnsu97WJBJ2WCZDAeest6Nq19kAHHSgV8ZJCXbJu4UIYPfrYr6unLuIdhbpk3aJFdpD0WBTqIt5R\nT12yav9+6NABtm2z1Rhrc+gQNGlit7XNjhFxmXrqEiivvgr9+h070MHWU2/SBD7/PHd1ieQrhbpk\nVdX89FR0sFTEGwp1yapFi45/kLSK+uoi3lCoS9bs3m1z04uLU++rUBfxhkJdsiYWgyFD4IQTUu+r\ns0pFvKFQl6xJt58OGqmLeEWhLlmTan56Ih0oFfGGQl2yYskSm6M+YEB6+2ukLuINhbp4rrISfvpT\nuOUWKChI72sU6iLeUKiL5xYuhA8/hEsvTf9rdKBUxBsKdfFUZSXcdBPcdhsUZnCxRPXURbyhUBdP\nPfOMreFSUpLZ16n9IuINLy88LY47cgRuvhl++UtokOFwQaEu4g2N1MUzM2dC8+Zw/vmZf6166iLe\n0NK74olDh6BPH5g+HUaMyPzrDxyAli3tVsvvilTLdOldhbrUW3k5TJpk959+uu7v07gxfPIJNG3q\nTV0i+UDrqUtOHTgAF1xg/fQnnqjfe6mvLlJ/CnWps337YPx4aNYMnnoqvYW7jkehLlJ/CnWpk927\nYexY6NQJHn/crl5UXzpYKlJ/6YT6w8DHwKrj7HM3sB5YCQz0oC4JqE8/hdJS6NkTzjgDHnkk/aUA\nUtFIXaT+0gn1R4Axx3n9PKAH0BO4Crjfg7okYLZuhWnToEcP2LIFli6F3/wm8/nox6OzSkXqL52T\nj14Gio7z+njg0fj914FWQEdsdC8hdPAgbN4Mr70GL74IL71kbZGSEnjzTSgqys731UhdpP68OKO0\nC7Al4fFWoCsK9UDbtAnWr7fwTt527YKuXWHQIBg+HH70I+jb19tReW0U6iL159UyAclzKGudkF5a\nWvrl/Wg0SjQa9ejbSyY++ghOOw2GDbNRd1GRHfTs1s3ud+rkXZ88E506wRtv5P77igRJLBYjFovV\n+evTndBeBMwF+tXy2gNADJgZf7wOGE7NkbpOPgqI9estxDds8LuSo61dC2PG2F8LOqtUxPhx8tEc\noGrl7GJgN2q9BFp5OTRq5HcVNfXubScxrV/vdyUi4ZVO+2UGNvJuh/XObwWqZiU/CMzDZsBsAPYC\nk70vU7x08GD9TxTKhkgEzjkHnn8eevXyuxqRcEon1Celsc/19S1EcieooQ5w9tnw5JNw3XV+VyIS\nTjqj1EHl5cEN9VGjIBazVR9FJHMKdQcdPBjMnjpAhw42C2fZMr8rEQknhbqDgtx+AWvBPP+831WI\nhJNC3UFhCPXnnvO7CpFwUqg7KMg9dYChQ2HVKtizx+9KRMJHoe6gIPfUAZo0gSFDYPFivysRCR+F\nuoOC3n4B9dVF6kqh7qCgt19AfXWRulKoOyjo7ReA/v2tp755s9+ViISLQt1BYWi/NGgAo0erBSOS\nKYW6g8IQ6qAWjEhdKNQdFIaeOtgyvAsXaskAkUwo1B0Uhp462EUzevSAl1/2uxKR8FCoOygs7ReA\nceNg7ly/qxAJD4W6g8LSfoHqUNdFs0TSo1B3UJhG6gMGwIED8M47flciEg4KdQeFpacOdjWkb38b\n/v53vysRCQeFuoPCNFIH9dVFMqFQd1CYeuoAI0fCihXwySd+VyISfAp1B4Wp/QK2amM0CvPn+12J\nSPAp1B0UtvYLqAUjki6FuoPCGOrnnw8LFujsUpFUFOoOCltPHaBzZzj1VFiyxO9KRIItnVAfA6wD\n1gNTa3m9HTAfeBtYDfzQq+IkO8LWU6+iFoxIaqlCvQC4Fwv2vsAkoE/SPtcDK4ABQBS4Cyj0tErx\nVBjbL6CzS0XSkSrUBwMbgM3AIWAmMCFpnw+BFvH7LYBdwGHvShSvhbH9AjBwoP1CWr3a70pEgitV\nqHcBtiQ83hp/LtF04KvAB8BK4D89q06yIqztl0gELrwQnnrK70pEgitVmySdP3RvwvrpUaA78Dzw\nNeDz5B1LS0u/vB+NRolGo+lVKZ4Ka/sFoKQEJk+Gn/3MQl4k38RiMWKxWJ2/PtXHohgoxXrqANOA\nCuDOhH3mAb8AquYlLMIOqC5Peq/KSjVDA6GwEPbvh4YN/a4kc5WV8JWvwLx5cPrpflcjkn0RG72k\nPYRJ1X5ZDvQEioBGwERgTtI+64DR8fsdgdOAsnQLkNw6cgQqKizYw0gtGJHjSxXqh7HZLQuANcCT\nwFpgSnwDuB34OtZPXwj8BNAqHQFVXm799DC3LkpKYNYsv6sQCaZcfrTVfgmA3butfbFnj9+V1F1F\nhf0bFiyAvn39rkYku7xuv0ieCet0xkQNGsB3v6sWjEhtFOqOCfPMl0QlJQp1kdoo1B0T1jnqyYYM\nsVbS2rV+VyISLAp1x+TLSF0tGJHaKdQdkw899SpqwYjUpFB3TL60XwC++U27xN2aNX5XIhIcCnXH\n5Ev7BawFc/HF8NhjflciEhwKdcfkU/sF4NJLLdSPHPG7EpFgUKg7Jp9G6mDrv3ToAIsX+12JSDAo\n1B2TTz31KpdeCn/6k99ViASDQt0x+TZSB5g0CebMgS++8LsSEf8p1B2Tbz11gI4dYehQmD3b70pE\n/KdQd0w+tl9ALRiRKgp1x+Rj+wVg/HhYsQK2bEm9r0g+U6g7Jh/bLwCNG9vFMx5/3O9KRPylUHdM\nvo7UoboFo2X7xWUKdcfka08dbNmAgwdhefLVcUUcolB3TD6P1CMRuPxymD7d70pE/KNQd0y+9tSr\nXHGFrdwY5sv1idSHQt0x+TxSBzjpJDj7bPjzn/2uRMQfCnXH5HNPvcrVV8MDD+iAqbhJoe6YfB+p\nA4wYYW2mV1/1uxKR3FOoOybfe+pgB0ynTLHRuohr0gn1McA6YD0w9Rj7RIEVwGog5kVhkh0utF8A\nLrsM5s6FnTv9rkQkt1KFegFwLxbsfYFJQJ+kfVoBvwPGAacDF3pco3jIhfYLQNu2MGEC/PGPflci\nklupQn0wsAHYDBwCZgITkvb5PvA0sDX+WGOjAHOh/VLl6qvhwQehosLvSkRyJ1WodwESl0jaGn8u\nUU+gDbAYWA5c4ll14jlXRuoAxcXQtCksWuR3JSK5U5ji9XQmhTUEzgBGAU2BpcBrWA/+KKWlpV/e\nj0ajRKPRNMsUr7jSUwc7YHrddXD33TZ3XSQMYrEYsViszl8fSfF6MVCK9dQBpgEVwJ0J+0wFmsT3\nA3gImA/MSnqvykpNHPZdv362kmH//n5Xkhv790NREbz4IvTu7Xc1IpmLRCKQOqu/lKr9shxrrxQB\njYCJwJykfZ4BhmIHVZsCZwJr0i1AcsulnjpAkyZwzTXwm9/4XYlIbqQK9cPA9cACLKifBNYCU+Ib\n2HTH+cA/gdeB6SjUA8ul9kuVa6+Fv/wFtm/3uxKR7Et7SO8BtV8CoFMnePNN6NzZ70py66qr7N+c\ncFhHJBQybb8o1B3Tti28+67dumTdOhg+HDZvtpaMSFh43VOXPOPSlMZEvXvD4MG6OLXkP43UHdOw\nIezd615fHWwGzFVXwdq10EDDGQkJjdTlmCoq4MgRC3YXDRsGzZvDs8/6XYlI9ijUHVJebiP0SC7/\nPguQSAT+53/gjju01rrkL4W6Q1ztpycqKYFdu7R0gOQvhbpDXJyjnqygAG6+2aY2arQu+Uih7hDX\nziY9losugh074IUX/K5ExHsKdYeo/WIKCuCWWzRal/ykUHeI2i/VLrrIlg3QaF3yjULdIRqpV1Nv\nXfKVQt0h6qkfrWq0vnix35WIeEeh7hCN1I9WWGij9Vtu0Whd8odC3SHqqdc0aRLs2QNzkq8SIBJS\nCnWHaKReU0EB/OpXMHUqHDrkdzUi9adQd4h66rUbMwa6doXp0/2uRKT+FOoO0Ui9dpEI/PrXcNtt\n8NlnflcjUj8KdYeop35sAwbAuefCnXem3lckyBTqDlH75fh+/nN44AHYutXvSkTqTqHuELVfju/k\nk2HKFPjf//W7EpG6U6g7RO2X1G68ERYsgGXL/K5EpG4U6g7RSD21Fi2sr37NNXaVKJGwUag7RD31\n9FxyiV327v77/a5EJHPphPoYYB2wHph6nP2+ARwGLvCgLskCjdTTE4nAfffBz34GH37odzUimUkV\n6gXAvViw9wUmAX2Osd+dwHwyuOq15JZ66unr2xeuvNKuaSoSJqlCfTCwAdgMHAJmAhNq2e8GYBaw\nw8vixFtqv2Tm5pthyRJdz1TCJVWodwG2JDzeGn8ueZ8JQFUHUuvdBZTaL5lp2hTuvhuuvdZ+diJh\nUJji9XQC+rfAjfF9Ixyn/VJaWvrl/Wg0SjQaTePtxStqv2Ru/Hj44x9tCYFf/MLvasQFsViMWCxW\n569P1f8uBkqxnjrANKAC659XKUt4n3bAPuDfgeTFTCsrtWi1r773Pfjud2HiRL8rCZePP4avfQ2e\neQbOPNPvasQ1kUgEMjhWmar9shzoCRQBjYCJ1AzrU4Fu8W0WcE0t+0gAqKdeNx07Whvmsstg/36/\nqxE5vlShfhi4HlgArAGeBNYCU+KbhIh66nX3ve/Zol9aQkCCLpfTD9V+8dnIkfDTn8KoUX5XEk67\ndkG/fjBzJgwb5nc14gqv2y+SR9R+qZ+2bW0Vx8mT4Ysv/K5GpHYKdYeo/VJ/48fbKP2GG/yuRKR2\nCnWHKNS9cc89sHQpPPqo35WI1KRQd4jmqXujWTN46ilbQmDNGr+rETmaQt0h6ql7p18/W6K3pAT2\n7vW7GpFqCnWHqP3ircmTYdAg9dclWBTqDlH7xVtVS/QuXWpLCYgEQaq1XySPaKTuvWbN4OmnIRqF\n3r2huNjvisR1Gqk7RD317OjbFx55xNbVef99v6sR1+mMUkdUVEBBgd1GdBmTrLjrLvjzn+GVV+DE\nE/2uRvJFpmeUKtQdceAAtGypdcGzqbISrrgCdu+GWbOggf4OFg9omQCplVov2ReJ2MWqd+zQwl/i\nH4W6I3SQNDdOOAFmz7aR+j33+F2NuEizXxyh6Yy50749PPccnHUWtGkDF1/sd0XiEoW6IzRSz62i\nIpg/35Y7btMGxo71uyJxhdovjlBPPfe++lX429/sikmvvup3NeIKhbojNFL3x5Ah8Nhj8J3vwBtv\n+F2NuECh7gj11P1z7rnw0ENw/vnw2mt+VyP5TqHuCLVf/DVunK2/Pn48LFnidzWSzxTqjlD7xX9j\nx9oZp9/5Drz4ot/VSL5SqDtCoR4M55xjF64uKYF58/yuRvKRQt0R6qkHx8iRMGcOXH45/OEPflcj\n+Ubz1B2hnnqwFBfDSy9ZS2bLFrj1Vi20Jt5Id6Q+BlgHrAem1vL6xcBK4J/AEqC/J9WJZ9R+CZ5e\nvWz++rPP2kJghw75XZHkg3RCvQC4Fwv2vsAkoE/SPmXAMCzM/w/4vYc1igfUfgmmjh1h8WLYvt36\n7Tt2+F2RhF06oT4Y2ABsBg4BM4EJSfssBfbE778OdPWoPvGI2i/B1awZPPOMnaj09a/Dm2/6XZGE\nWTqh3gXYkvB4a/y5Y7kC0HH9gFH7JdgKCuD22+1CG2PGwJ/+5HdFElbpHCjN5MoWI4DLgW/V9mJp\naemX96PRKNFoNIO3lvpQqIfDhRdCnz7wb/9mZ5/edRc0aeJ3VZJLsViMWCxW569P53h7MVCK9dQB\npgEVwJ1J+/UHZsf321DL++jKRz6qml2R8HtVAmz3brj6avjXv2DGDDj9dL8rEr9k48pHy4GeQBHQ\nCJgIzEna5xQs0H9A7YEuPlNPPVxatbIw/+//hhEj4L777HJ5IqmkE+qHgeuBBcAa4ElgLTAlvgHc\nArQG7gdWAMs8r1TqRe2X8IlE4Ic/tLViHn7Y1o3Zts3vqiTo0p2n/g/gNKAH8Mv4cw/GN4ArgbbA\nwPg22MMaxQMK9fCqms8+aBAMGGArPmrULseiZQIcUV6ueeph1qiRHQ9ZtAgeeADOPhvKyvyuSoJI\noe4IjdTzQ//+Nivm3HPhG9+A226D/fv9rkqCRKHuCIV6/igshB//2E5SWrmy+rJ5askIKNSdoWUC\n8k9RETz9NPz+93DTTTZ6X7nS76rEbwp1R2hKY/4aPdrCfNw4C/aLL4aNG/2uSvyiUHeE2i/5rWFD\nuOEGWL8eeveGM8+E667TFEgXKdQdoVB3Q/PmcPPNsG6dLS/Qrx9MmaKZMi5RqDtCPXW3tGsHv/41\nvPMOtG8PgwfDJZfAqlV+VybZplB3hHrqbmrfHn7+c+ux9+1rK0COHGlL/R454nd1kg0KdUeo/eK2\nli1h2jTYtAmuvNKW+e3Vy0bzujBHflGoO0LtFwH7f+D734fXX4fHH4fVq6FnTygpgQULNHrPB7m8\n1K2W3vVRly6wbJndiiTas8dWhHzoIbus3qRJNi2yXz9dDDsIMl16V6HuiHbtbEZEu3Z+VyJBtmoV\nPPGEbc2b26j+wgutVSP+UKhLrZo3hw8+sFuRVCoqbGXIGTPgr3+F1q3hggtsGzBAI/hcUqhLrRo1\ngs8/18FSyVxFhbXuZs+27cABGDvWttGjoUULvyvMbwp1qaGyEho0sA+nRlhSH5WVdtbqP/4B8+bZ\naH7gQJsmOWqUncmqA/LeUqhLDQcP2mjq4EG/K5F8s3cvvPIKvPCCrfX+7rtQXAxDh8JZZ9lJTyee\n6HeV4aZQlxo++wy6drVbkWz65BO7/N4rr9i2ciX06WMj+MGDbevVy/5ylPQo1KWGnTttkaedO/2u\nRFyzf7+t+/7GGzY3ftky2LXLDrYOHFi99e6tts2xKNSlhm3bbISkFfskCHbuhBUr4O237XbFCti8\nGbp3h9NPt61vXwv6Hj0U9gp1qaGszGYpaKU+CaoDB+w8itWrba78unW2vfcenHKKtWx69Kjeune3\n512YzZVpqBdmrxQJCi0RIEHXuLG1ZAYMOPr58nJbjGz9etiwAdasgTlz7LmtW23Bsm7d7CpQJ59s\nQX/yybZ16QJt2rg34yudUB8D/BYoAB4C7qxln7uBscA+4IfACo/qEw9ohUYJq0aN7EBrnz41Xzt8\n2FqKmzfbtmULvPWWrUD5/vt2st3+/dC5s22dOsFJJ9nWqRN06GBb+/Z2my+zdFKFegFwLzAa2Aa8\nAcwB1ibscx7QA+gJnAncDxR7XmmOxGIxotGo32WklEmdfq7QGIafZxhqBNWZrLAQvvIV24YPr32f\nffss3Ldtg48+qt5eegnWro1RURFl+/bqlSrbtrWlNNq1s/utW9tov00bu9+qVfXWsmX11rBh1v+5\naUsV6oOBDcDm+OOZwASODvXxwKPx+68DrYCOwMeeVZlD+fjBUagfXxhqBNVZF02bVvfhk5WWxigt\njQJ2UtW+fTYzZ+fO6ttPP7Vpmh98YP3+3buP3j77zBZEa9jQzgVp0cKW4kjcmjWzvwKqbk880eqq\num3a1K5SVXVbtTVubLeZShXqXYAtCY+3YqPxVPt0JaShno/Ky9VTFzmeSKQ6cE85JbOvray0Ns+e\nPbYUR+L2xRd2gtYXX9i2Zw98+KE9V7Xt32/bvn22HThgW9XzmUoV6ulOV0k+FFHr140bl+a7+eid\nd2xebdBlUuf27VqdUSRbIpHqEXenTtl5/4z2T/F6MVCKHSwFmAZUcPTB0geAGNaaAVgHDKfmSH0D\n0D2z8kREnLcRO27picL4GxYBjYC3geTj0OcB8+L3i4HXvPrmIiLivbHAO9hIe1r8uSnxrcq98ddX\nAmfktDoREREREambMViffT0w1edaEj2M9f1XJTzXBngeeBd4Dpue6beTgcXAv4DVwH/Enw9arY2x\nKa1vA2uAX8afD1qdVQqwk+Tmxh8Hsc7NwD+xOpfFnwtana2AWdg05zXY7Lig1Xga9jOs2vZgn6Og\n1QnWDfkXlktPACcQsDoLsLZMEdCQ2nvyfjkLGMjRof4r4Cfx+1OBO3JdVC1OAqpOnm6GtcL6EMxa\nm8ZvC7FjK0MJZp0APwIex06mg2DWuQn7QCcKWp2PApfH7xcCLQlejYkaAB9ig6Wg1VkElGFBDvAk\ncBkBq3MIMD/h8Y3xLSiKODrU12EnToGF6bpcF5SGv2Fn+Aa51qbY2cdfJZh1dgUWAiOoHqkHsc5N\nQNuk54JUZ0sshJIFqcZk5wAvx+8Hrc422KCtNfYLci5wNgGr80JgesLjHwD3+FRLbYo4OtQ/Tbgf\nSXocBEXAe0BzgllrA+yvsc+x0QUEs86nsL/ShlMd6kGsswxrFywH/j3+XJDqHIC13B4B3sI+6ycS\nrBqTPQxcG78fxDqvwj4/24HH4s9lVGe2rz8S5rV2KwlW/c2Ap4H/xP6jJwpKrRXYB70rMAwbCScK\nQp3fxj4wKzj2eRpBqBPgW9gvn7HAdVjLMJHfdRZis93ui9/upeZf4n7XmKgRMA77pZ4sCHV2B/4L\nG7x1xj7zP0jaJ2Wd2Q71bVjvqsrJ2DICQfUx9ucNQCfswx8EDbFAfwxrv0BwawU7EPUsMIjg1flN\nbL2iTcAMYCT2cw1anWC9X4AdwF+xtZiCVOfW+PZG/PEsLNw/Ijg1JhoLvIn9PCFYP0uArwOvAruA\nw8BsrIWd0c8z26G+HFu9sQj7LTmR6gNTQTQHOzBB/PZvx9k3VyLAH7CZBb9NeD5otbaj+qh8E6wX\nuILg1XkTNrjoBlwEvABcQvDqbIq12cBaGudgrcIg1fkRtu5Tr/jj0djMjbkEp8ZEk7Bf5FWC9LME\n65UXY5+fCPbzXEMAf561nbwUBDOAD4By7H/MydiBioUEZOpQ3FCsrfE21VOyxhC8WvthfdW3sWl4\nP44/H7Q6Ew2nepARtDq7YT/Lt7GprFWfnaDV+TVspL4SG1m2JHg1gv1i3En1L0oIZp0/oXpK46PY\nX+lBrFNEREREREREREREREREREREREREREREREREwub/AU1lgzCZ356cAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "r_ir = relativistic_ir(dt=0.125)\n", "plt.plot(r_ir)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }